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1. INTRODUCTION TO THE NOTES 


Quantum information is an emerging field which has attracted a lot of attention in the last fifteen years or so. It is a 
broad subject which covers from the most applied questions (e.g., how to build quantum computers or secure crypto¬ 
graphic systems), to the most theoretical problems concerning the formalism of quantum mechanics, its complexity, 
and its potential to go beyond classical physics. 

I have written these notes as an introduction to the formalism of quantum information, trying to show how it is a 
natural progression from an undergraduate course on quantum mechanics (just as I felt it was for me). Hence, these 
are mainly notes about the language of quantum information and the meaning of its terms, not about its applications. 

Not showing many examples of the application of the tools that are introduced in the notes, the reader might 
feel that some definitions are just ‘highly complicated ways of writing simple quantum mechanical expressions’, even 
finding hard to understand the motivation behind introducing concepts like ’entanglement’ or ‘quantum operations’; 
questions along the line of “so what?” might cross the reader’s mind, and under such circumstances I can only ask 
for his/her patience: this is a work in progress and I hope to keep improving the notes with practical examples, 
applications, and implementations... but for the time being, just trust me when I say that there is a proper reason for 
introducing every concept appearing in the notes, as well as for writing each expression the way it is written. 

Also, let me remark that these notes do not intend to be a comprehensive and exhaustive review about the field, 
but just a logical sequence of ideas, concepts and definitions, building up towards the topics I’ve worked on within the 
continuous variables formalism, which actually takes over most of the notes. In this way, most of the results are stated 
without proof, and I haven’t been careful with references; as I said, this is a work in progress, and I intend to add 
such things in the future. In any case, most of the topics that I introduce have been thoroughly covered in previous 
books and reviews mil El HE], and with this notes I only intend to give my point of view, as well as to introduce the 
different topics and objects in exactly the way I work with and think about them, hoping it will be useful, especially 
for the people who are starting in the field. 

This notes originated in 2010 as a way of studying for my close collaboration with Raiil Garcia-Patron and Nicolas 
J. Cerf, to whom I am deeply grateful for guidance and incredibly clarifying discussions. They benefited greatly as 
well from the weeks I spent in Doha with Hyunchul Nha and his group in 2012, especially from discussions with him 
and Ho-Joon Kim. 

Munich, April 21, 2015. 


Carlos Navarrete-Benlloch 
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1. Introduction to the notes 



2. QUANTUM MECHANICS IN AN ISOLATED SYSTEM 


The intention of this initial chapter is to introduce the axioms of quantum mechanics as we will use them during 
these notes. These are introduced from a very pragmatic viewpoint, not trying to justify them or motivate them 
from a physical perspective, for a reasoned presentation of them I recommend reading [B] (Appendix A, as well as the 
references therein), as the tone and notation of these notes are fully correlated with that reference. The mathematics 
behind quantum mechanics (Hilbert spaces) are assumed to be known, although I provide a summary of them in 
Appendix m as a reminder, and in order to fix the notation. 

2.1 The state of the system 

In the mathematical framework of quantum mechanics, a Hilbert space H is associated to any physical system. The 
state of the system is completely specified by a density operator p, that is, a self-adjoint (p = pi), positive (('0|p|'0) ^ 0 
V|'0) G %), and unit-trace (tr{p} = 1) operator. 

Density operators can always be written as a convex mixture of projectors, that is, 

N 

P = ( 2 . 1 ) 

k^l 

where {wk}k=i, 2 ,...,N is a probability distribution {wk > 0 V/c and = 1) ^'^*1 states {\‘Pk)}k=i, 2 ,...,N are 

normalized, but may not be orthogonal. The set {wk, \'Pk)}k=i, 2 ,...,N is known as an ensemble deeomposition of the 
state p, and it is not unique (in the sense that different ensembles can lead to the same p). Not that N does not have 
to coincide with the dimension of the Hilbert space H. 

It can be proved that two ensembles {wk^\^Pk)'\k=l ,2 ...,n and {t'i, |'0i)}j=i,2,...,M give rise to the same density 
operator p if and only if there exists a left-unitary matri}|3 Li such that [T] 

y/^\^k) = ( 2 . 2 ) 

i 

where ii k f then |fc — /| zeros must be included in the ensemble with less states, so that lA is & square matrix. 

When only one state |<p) contributes to the mixture, then p = |<p)(<p| and we say that the state is pure] otherwise, 
the state is mixed. A necessary and sufficient condition for p to be pure is jT = p. 

In the next chapters we will learn that the mixedness of a state always comes from the fact that some of the 
information of the system has been lost to some other inaccessible system with which it has interacted for a while 
before becoming isolated. In other words, the state of a system is pure only when it has no correlations at all with 
other systems. 


2.2 Evolution of the system 

Quantum mechanics contemplates two different types of evolution of a closed system which is initially in some state 

/ 3 : 

• Unitary. If the system is allowed to evolve during a time t without ‘observing it’, its state changes to 

p' = t]pij\ (2.3) 

where 

U = eyip{Ht/ih), (2.4) 

is a unitary operator, being H a self-adjoint operator corresponding to the Hamiltonian of the system. 

^ W is left-unitary if lA^lA = X but tUA^ ^ X, where X is the identity matrix of the corresponding dimension. It is easy to prove that 
finite-dimensional left-unitary matrices are unitary. 
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2. Quantum mechanics in an isolated system 


• Projective measurements. When one measures an observable with associated self-adjoint operator Y = 
yjPj^ where Pj = \yj){yj \ are the projectors onto its eigenvectors^ and yj G M, the outcome yj appears with 
probability 

Pj = tr{/5Pj} = (yjlplyj)- (2.5) 

If the outcome is recorded by the observer, the state of the system changes (collapses) to 

P'j = Pj^PoPP] = \V]){y]V ( 2 - 6 ) 

and the measurement is called selective. If, on the other hand, the outcome of the measurement is unknown to 
the observer, the system is left in the mixture 

P' = Y1 Po 1%) (yj I = X! P^iPPj ’ (2.7) 

3 3 

and the measurement is called non-selective. 

It might seem strange to see that the state of the system is different after the measurement depending on whether 
the observer learns the outcome or not; however, note that the state (1^ is statistically indistinguishable from 
keeping track of the outcomes at any idividual realization of the measurement, and assume that the state 
collapses after each measurement to the corresponding state (ESI). That is, on what concerns to the statistics of 
any experiment, selective and non-selective measurements are equivalent. 

Note that these two types of evolution are fundamentally different: while unitary evolution is reversible, projective 
measurements are not. 


2.3 The von Neumann entropy 


The von Neumann entropy of a state p is defined as 

5[p] =-tr{/51ogp}. (2.8) 

Given the diagonal representation of the state 

d 

P = '^>^3p3){r3\^ (2-9) 

f=l 

where the set {\rj)}j-i^ 2 ,...,d forms an orthonormal basis of the Hilbert ^ace "H (which therefore has dimension d), 
the von Neumann entropy is just the Shannon entropy of the distributioijj A = col(Ai, A 2 ,..., Xd), that is. 




d 

^AjlogAj. 

j=i 


( 2 . 10 ) 


For a pure state |^)(^| the entropy is zero, while it has a miximum \ogd for the maximally mixed state 



d 


( 2 . 11 ) 


Hence, the entropy can be understood as a measure of the mixedness of the statcQ. Note that when the state of the 
system is /3mm and a given observable Y acting on H is measured, all its eigenvalues are equally likely to appear as 
an outcome of the measurement, that is. 


Pj = {yj\p\yj) = 1/d Vj. [flat distribution] (2.12) 

This reinforces the interpretation of mixedness as due to some kind of information loss. 

The entropy is a concave functional of density operators, that is, given the convex mixture of states 

P = '^WkPk: (2-13) 

k 

^ We will only consider the case of observables having non-degenerate eigenvalues. 

^ In the following we use the notations p = co\{pi,p2<, •■■■,Pd) and for a probability distribution interchangeably. 

The case d oo will be discussed later when studying in detail infinite-dimensional spaces. 




2.3. The von Neumann entropy 
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where {wk}k=i.2,...,N is a probability distribution and {pk}k=i,2,...,N is a set of density operators, we have 

N 

sM > E WkS[pk\. (2.14) 

k^l 

This is easily proved by writing each density matrix in its diagonal form, and using the concavity of the —xloga; 
function. 

Finally, it is interesting to note that the entropy does not change with unitary evolution, can only increase with 
non-selective projective measurements, and can only decrease with selective projective measurements (indeed, when 
no degeneracies are present, it collapses to zero, as the state becomes pure). 



6 


2. Quantum mechanics in an isolated system 



3. BIPARTITE SYSTEMS AND ENTANGLEMENT 


3.1 Entangled states 

Consider two systems A and B (named after Alice and Bob, two observers which are able to interact locally with their 
respective system) with associated Hilbert spaces Ha and Hb- Imagine that the systems A and B interact during 
some time in such a way that they cannot be described anymore by independent states pA and pB acting on Ha and 
Hb, respectively, but by a state pab acting on the joint space Ha®Hb- The question is, is it possible to reproduce the 
statistics of the measurements performed by Alice on system A via some state pA defined in Ha only? This question 
has a positive and unique answer: this state is given by the reduced density operator 

PA = Hb{pab}, (3.1) 

that is, by performing the partial trac^il over system’s B subspace onto the joint state. Of course, the same applies to 
Bob, whose system is described locally by the reduced density operator pB = HaIpab}- 

If the state of the joint system is of the type = Pa® Pb, that is, a tensor product of two arbitrary density 
operators, the actions performed by Alice on system A won’t affect Bob’s system, the statistics of which are given by 
Pb, no matter the actual state pA- In this case A and B are uncorrelated. For any other type of joint state, A and B 
will share some kind of correlation. 

Correlations are not strange in classical mechanics; hence, a problem of paramount relevance in quantum mechanics 
is to understand which type of correlations can appear at a classical level, and which are purely quantum. This is 
because only if the latter are present, one can expect to hnd quantum mechanical effects such as violation of local 
realism, computational algorithms exponentially faster than the fastest ones known in classical computer science, or 
unconditionally secure cryptography. 

Intuitively, the state of the system will induce only classical correlations between A and B if and only if it can be 
written as 



fc 


where the p^^’^s are density operators and {wk}k=i, 2 ,... is a probability distribution, because it means that Alice and 
Bob can share a classical machine which randomly picks a value of k according to the distribution {wk}k to trigger 

'^(k) '^(k) 

the preparation of the states p\' and pg, what can be done locally, and hence cannot induce further correlations. If 
the process is automatized so that Alice and Bob do not learn the outcome k of the random number generator, the 
mixture p\g is prepared. In other words, the state does not contain quantum correlations if it can be prepared using 
only local operations and classical communicatioiH. 

There is yet another way of justifying that states which cannot be written in the separable form (ESI) will make 
A and B share quantum correlations. The idea is that one of the fundamental differences between classical and 
quantum mechanics is the superposition principle, that is, interferences. Hence, it is intuitive that correlations will 
have a quantum nature when the joint state of the systems exploits the concept of “superposition of joint states”, that 
is, when it cannot be written as a tensor product of two independent states of the systems, or as a purely classical 
statistical mixture of these, which corresponds exactly to dSED. 

States of the type p^^ are called separable. Any non-separable state will induce quantum correlations between A 
and B] these correlations which cannot be generated by classical means are known as entanglement, and states which 
are not separable are called entangled states. 

^ Given an orthonormal basis {|bj)}j of 'Hb, this is defined by 

tr{PAfl} = ^(fedPAsIfej), (3.2) 

j 


which is indeed an operator acting on I-La- 

^ We will learn the exact meaning of this class of operations in the next chapter. 



3. Bipartite systems and entanglement 


3.2 Characterizing and quantifying entanglement 

In general, given a mixed state pab acting on Ha ®'Hb, it is hard to find out whether it is separable or not, the 
difficulty coming from distinguishing between quantum and classical correlations. Indeed, the best known criterion for 
separability, the Peres-Horodecki criterion, yields only necessary and sufficient conditions when dAXdp < G (note that 
this includes the case of two qubits), and for a reduced class of states in infinite-dimensional Hilbert spaces (Gaussian 
states, see Chapter [6l). This criterion states that a necessary condition for the separability of a density operator is 
that it remains being positive after the operation of partial transposition, that is, given 

Pab — ^ ( Pil,km\t^j (^■^) 

jklm 

where {\aj,hi) = |aj) 0 \bi)}j,i is an orthonormal basis of Ha 

Pab — ^ ^ I 5 1 1 5 (^■^) 

jklm 

is a positive operator. We will learn more about this criterion and some more when studying infinite dimensional 
Hilbert spaces in Chapter [51 

A very different problem is that of quantifying the level of correlations present in the state, and more importantly, 
how much of these correspond to entanglement. Even though we pretty much understand the conditions that a 
proper entanglement measure E[pab] must satisfy, we have not found a completely satisfactory one for general states 
[3] (either they do not satisfy all the conditions, or/and can only be efficiently computed for restricted classes of 
states). It is not the intention of these introductory notes to explain all these measures and up to what point they are 
satisfactory, but let me spend a few words on this topic for completeness (see [2112] for more details). 

The basic conditions that a good entanglement measure E[pab] should satisfy are quite intuitive: 

1. E[pab] is positive definite and equal to zero for separable states. 

2. Given the mixture pab = where {pj}j is a probability distribution and {pj}j are density operators 

acting on Ha ® Hb, E[pab] < 'Bj which is as to say that the entanglement of a collection of states 

cannot be increased by not knowing which of them has been prepared. 

3. At least on average, the entanglement can only decrease when Alice and Bob apply protocols involving only local 
operations and classical communication. 

These three conditions define what is known as an entanglement monotone. By themselves, they are not enough 
to define a unique entanglement measure even for pure states. However, by adding two more conditions known as 
weak additivity and weak continuity |3], which find an intuitive justification in the asymptotic limit of having infinitely 
many copies of the state, it is possible to prove that the entanglement entropy, which measures how mixed is left 
the reduced density operator of one of the parties after tracing out the other, is the unique entanglement measure of 
pure states pab = |V')ab(' 0I) given the reduced density operator pA = trs{|'0)AB('0|} or ps = tr^{|'0)^B('0|}, this 
entanglement measure can be evaluated as 

E[\iIj)ab] = S[pa] = S[pb], (3.6) 

where the equality of the von Neumann entropies of the reduced states will be clear after the following section. Hence, 
the problem of quantifying the entanglement is basically solved for pure bipartite states. The pure states whose 
corresponding reduced states are maximally mixed are known as maximally entangled states. 

In the case of mixed states, the entanglement entropy is not even an entanglement monotone, as the entanglement 
entropy of the separable state pa® Pb is just the entropy of the mixed states pA or pp, which is not zero except for 
pure states (note that, furthermore, it depends on which mode is traced out). However, there are many quantifiers 
which are entanglement monotones: 

• Possibly the most natural entanglement measure for mixed states is the distillable entanglement. Suppose that 
we give N copies of the mixed state pab to Alice and Bob; the process of distillation refers to the conversion of 
these copies to copies of maximally entangled states via protocols involving only local operations and classical 
communication. The distillable entanglement is defined as the maximum number of maximally entangled states 
that can be distilled from infinitely many copies of pab- Apart from an entanglement monotone, it can also be 
shown to satisfy the weak additivity and weak continuity conditions, and to be equal to the entanglement entropy 
for pure states. Its drawback is that it requires a maximization over all the possible distillation protocols, and 
it is therefore very hard to evaluate (paraphrasing |3]: it is a problem ranging from “difficult to hopeless”.) 



3.3. Schmidt decomposition and purifications 
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• The entanglement of formation can be seen as the dual of the distillable entanglement: it measures the number 
of maximally entangled states that are needed to prepare infinitely many copies of the mixed state. It can be 
evaluated as 

Ef[pab]= min ^ Wfc£;[|V’fe)], (3.7) 

{wk,\^k)}k ^ 

where the minimization is performed over all the possible ensemble decompositions of pab-, what makes the 
measure very hard to evaluate again. Nevertheless, closed formulas have been obtained for the entanglement of 
formation of the general state of two qubits x ds = 2 x 2), as well as for reduced classes of higher-dimensional 
bipartite states with strong symmetries. 

• There is an entanglement monotone which can be evaluated fairly efficiently, as it does not require any opti¬ 
mization procedure: the logarithmic negativity. In loose terms, it quantifies how much the state pab violates 
the Peres-Horodecki criterion via 


En[pab] = log||p5i3lli 


log 


1+I “ 

3 


(3.8) 


where {Aj}j are the eigenvalues of p^g, and ||A||i = trvXIt denotes the so-called trace norm. The problem 
with this measure is that it does not collapse to the entanglement entropy for pure states, and is not weakly 
additive in general. 


3.3 Schmidt decomposition and purifications 

It is possible to show [T] that any pure bipartite state can be written in the form 

d 

|^)ab = ^ ® (3-9) 

3 = 1 

where {Aj}j=i_ 2 ,...,d is a probability distribution, and {|wj)}j=i, 2 ,...,d and {\vj)fj-\^ 2 ,...,d are orthonormal bases of T-La 
and %B , respectively, which we assume to have the same dimension without loss of generality. With the state written 
in this form (known as Schmidt decomposition), it is completely trivial to evaluate the entanglement entropy: as the 
reduced density operators are diagonal, that is, 

d d 

Pa=J 2 ^3\^3){uj\, and PB I’ (3-10) 

1=1 i=i 

their von Neumann entropies give 

d 

E[\d)AB] = 

3 = 1 

The Schmidt decomposition allows us to introduce the concept of purification: given a system A with associated 
Hilbert space "Ha in a mixed state with diagonal representation pA = can always introduce another 

system B laying in a Hilbert space T-Lb with the same dimension as T-La, and with an orthonormal basis {Ir’j)}^, and 
understand the mixed state pA as a reduction of the pure entangled state |V')ab = 'Yfj ® bi)- Hence, a 

classical mixture of states can always be transformed into a pure entangled state in a doubled Hilbert space. 
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3. Bipartite systems and entanglement 



4. QUANTUM OPERATIONS 


4.1 Quantum operations as evolution of a reduced system 


Consider a system Q with associated Hilbert space Hq subject to the actions of an experimentalist named Quinn; 
quantum operations appear as an answer to the following question: are unitary transformations and projective mea¬ 
surements the only type of operations that Quinn is allowed to perform onto the system Q? 

The answer to this question is ‘no’: Quinn can always append an auxiliary system L with associated Hilbert space 
Hl, apply unitaries and projective measurements on the joint system, and finally dismiss (trace out) the appended 
system L. Quantum operations correspond to these reduced unitaries and projective measurements as felt by the 
system Q alone. 

As we are going to prove in the reminding of the section [T] , any quantum operation can be represented by a map 
of the type 

K 





(4.1) 

where K can be selected at will by Quinn, 

and where the only restriction on the operators {Ek\k=i, 2 ,. 

,,^K (which act 

onto Hq) is that they must satisfy 

( 

■ 

\ 




iZ^lEk 

k—1 

/5 i < 1 V/5, 

(4.2) 

what is denoted symbolically by 


ic 





Y,ElEk<i-, 

(4.3) 


k=l 


in other words, the map can be trace preserving or trace decreasing. If the state of the system Q was pQ prior to the 
quantum operation, it becomes 


-/ ^ ^[pq] 

^ tr{£:[pQ]}’ 


(4.4) 


after it. The operators {Ek}k=i,2,...,K are known as Kraus operators, and the expression (SH) as an operator-sum 
representation of the quantum operation. 


Let’s see how by applying unitaries and measurements onto the joint Hilbert space T-Lq®T-Ll, the reduced dynamics 
of system Q is described by a map of the type 611). Suppose that the initial state of the joint system is 

PQL = PQ®\yL){yL\] (4.5) 

note that assuming that system L is in a pure state is still completely general, because if it is in a mixed state, this 
can be purified by introducing another extra system, and renaming the joint system as L. 

If Quinn applies a global unitary U, the state evolves to 

p'QL = U{pQ®\yL){yL\)U'^ ■ (4.6) 

Introducing now an orthonormal basis in given by , the reduced state of the system Q is written as 

di, di, 

p'q = ^^lVp'ql) = Y.^lk\U{PQ ® WL){yL\)U^\lk) =Y,EkPQEl (4.7) 

fc=l k=l 

where the operators Ek = {IklU^L) act onto Hq. Note that trQ{;OQ} = 1, and hence, the reduced unitary corresponds 
to a trace-preserving map like 6ID with a number of Kraus operators K given by the dimension of the Hilbert space 
of the appended system L. 
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4. Quantum operations 


If, on the other hand, Quinn applies a non-selective measurement of a joint observable with diagonal representation 


dq dL 

Y=J2yjPj, ( 4 . 8 ) 

i=i 

where {Pj = \yj){yj\}j=i, 2 ,...,dQ-dL ^^re the projectors onto its eigenvectors {\yj) € Hq 'Si 'HL}j=i, 2 ,...,dQ-dL^ state 
evolves to 

dq-dL 

PQL= P3ih®\^L)i^L\)Pj- (4.9) 

i=i 

The reduced state of system Q can be written then as 


PQ 


U.Q-U.L aL U-QU-L aL 

= ^^l{p'ql}= Y. Y^h\P,{pQ®WL){'PL\)Pj\lk) = Y Y^k''PQE\^'’^ 

j—1 k—1 j—1 k—1 


( 4 . 10 ) 


where = {lk\Pj\'PL)■ Again, it is immediate to check that tcglpg} = 1, and hence, the reduced non-selective 
measurement is a trace-preserving quantum operation with a number of Kraus operators given by dg ■ <Pj^. 

The reduced dynamics of a selective measurement of the joint observable K is a little more subtle. Assume that 
after the measurement Quinn obtains the outcome so that, accordingly, the state pql collapses to 


Pql=Pj ^PjiPQ ® Wl){'Pl\)P], 


(4.11) 


where 

P3 = tr{PjPQL}, (4.12) 

is the probability for the outcome yj to appear after the measurement. In this case it is easy to rewrite the reduced 
state of system Q as 

Y ( 4 . 13 ) 

k^l 

and the corresponding probability as 


P3 = IrQ 




.fc=i 



< 1, 


(4.14) 


where = {lk\Pj\'PL)- Hence, similarly to the previous cases, the reduced selective measurement is described by a 
map of the type 611), but the map is trace decreasing in general. 


Indeed, an important theorem states that any trace preserving quantum operation can be written as the reduced 
evolution of system Q after the application of a joint unitary transformation acting onto the Hilbert space T-Lq SPl- 
On the other hand, trace decreasing operations require an extra measurement of a joint observable, and the desired 
quantum operation is accomplished only when a particular outcome appears in the measurement, what happens 
with a probability equal to the trace of the map (it requires then post-selection). Hence, while trace preserving 
operations can be implemented deterministically—that is, at will by ‘pushing a button’—trace decreasing operations 
can only be applied probabilistically—meaning that we need to wait for the appareance of a particular outcome in 
some measurement—. We will sometimes refer to system L as the environment, and a representation of the quantum 
operation in terms of reduced unitaries and measurements acting onto the system and the environment is known as a 
Stinespring dilation. 

Finally, it is possible to show [T] that two sets of Kraus operators {Ek}k=i,2,...,K and {Pi}i=i,2 . l lead to the 

same quantum operation if there exists a left-unitary matrix U for which 

Ek=Y^kiPu (4.15) 

I 

where if K ^ L, then \K — L\ zeros must be included in the set with less Kraus operators, so that ZY is a square matrix. 





4.2. Generalized measurements and POVMs 
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4.2 Generalized measurements and POVMs 

From the previous discussion it should be clear that the most general measurement that one can perform on a system is 
completely described by a collection of trace-decreasing quantum operations {£j}j=i, 2 ,...,.j>i, which forms a complete 
set, that is, 


'^tT{£j[p]} = iyp. (4.16) 

The generalized measurements whose associated set of quantum operations {£j}j 
single Kraus operator, that is, 

£j[p] = Ej0. 

with 

,7 

tT{EjEjp} < 1 V/3 and ~ G 

7 = 1 

are very special, because it can be shown that its simplest Stinespring dilation does not require a joint unitary, just a 

joint projective measurement. The set {flj = EjEj}j^i^ 2 , _ r is known as a positive operator-valued measure (POVM), 

while the operators are called the measurement operators. 

This generalized measurements are the closest ones to projective measurements; the POVM {IIj}j=i, 2 ....,j plays 
the role of the spectral decomposition of the measured observable, while the measurement operators {Ej}j-i^ 2 ,...,j 
(uniquely dehned from the POVM up to a left-unitary transformation) play the role of the projectors. The probability 
of observing the outcome ‘j’ is 

Pj = (4-19) 

if the system was in the state p prior to the measurement, after which it collapses to the state 

p,=p-^E,pE}, (4.20) 

if the measurement is selective, or the state 

J 

p' = J2e,0^, (4.21) 

7 = 1 

if it is non-selective. 

As a simple application of POVM-based measurements, consider the following problem. Suppose that someone 
picks one state out of the set {\'Pi)dP‘ 2 )} and asks us to hnd with a single measurement which one did he/she 
picked. If the states are orthogonal, this is trivial: we make a projective measurement defined by the projectors 
{Pi = \pi){ipi\, P 2 = \‘P 2 ){'p 2 W, and check which outcome appeared. The problem is that it is simple to prove that 
when the states are not orthogonal, there is no strategy based on projective measurements allowing us to determine 
which state was given to us. However, we can design a strategy based on POVMs which will allow us to perform the 
needed task, although it does not work all the time. 

Consider the POVM {Hi = ki{i — \ipi){ipi\),ti 2 = k 2 {i — |v72)(</^2|),Ha = 7 — Hi — { 12 }. Suppose that we get the 
outcome ‘1’; then, we know for sure that we got the state \<p 2 ) as the probability of observing ‘1’ when the state is 
\ipi) is zero, that is, (:pi|lli= 0. The opposite happens when we get the outcome ‘2’, we know for sure that \Lpi) 
was given to us, because {p 2 \it. 2 \^P 2 ) = 0. Finally, when we get the outcome ‘3’ we don’t know which state we had, 
but at least we never make a misidentification of the state. 

4.3 Local operations and classical communication protocols 

Finally in this chapter on quantum operations, we consider a very important class of operations performed onto a 
bipartite system as the one considered in Chapter 2. 

Suppose that Alice and Bob are in distant locations, so that one does not have access to the part of the system 
belonging to the other. In this scenario, it is natural to think that the most general class of operations that can be 
performed on the joint system are local operations (arbitrary operations acting only on A or i? independently) in a 
correlated fashion (Alice and Bob can communicate by phone to decide together what to do). Quantum operations 
of this kind are known as loeal operations and classieal communication (LOCC) protocols, and play a central role in 
many problems of quantum information (we already saw one, the characterization of entangled states). 


=i, 2 ....,j>i are described by a 

(4.17) 

(4.18) 
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4. Quantum operations 


As an example, consider the following prototypical protocol: Bob performs a generalized measurement described by 
the set of trace decreasing quantum operations, and communicates the outcome, say ‘j’, to Alice, who 

applies a trace-preserving quantum operation from a pre-agreed set in one-to-one correspondence 

with the possible outcomes of Bob’s measurements. Let us denote by {Bl^^}k=i, 2 ,...,Kj the Kraus operators associated 
to £g'^ (note that when Kj = 1 Vj, Bob’s measurement is a POVM-based measurement), and by {Am}m=i,2,...,Mj 
the ones associated to £^f^. The resulting possible maps will be given by 

Kj Mj 

4 occ = £a [£b^[ p]] = EE B^k^ypii® ® y- (4.22) 

k—1 m—1 

We will find this type of LOCC protocols many times along the notes, which we will call one-way-direct LOCC 
protocols. 

Note that trace preserving LOCC protocols can only decrease the entanglement of the state shared by Alice and 
Bob, as intuition says. What is a little more surprising is that trace-decreasing LOCC protocols can enhance the 
entanglement (we shall find one example of this when studying photon addition and subtraction), what is a further 
example of how much counter-intuitive quantum mechanics can be. Of course, a complete set of trace decreasing 
LOCC maps can only decrease the entanglement on average, showing that any local operation able to enhance the 
entanglement is intrinsically probabilistic. 



5. MAJORIZATION IN QUANTUM MECHANICS 


In this section we introduce a relation between probability distributions called majorization [7], and connect it to a 
couple of important questions in quantum mechanics [T]. 

5.1 The concept of majorization 

Majorization appeared as a way to order vectors in terms of their disorder, in an effort to understand when one 
probability distribution can be built from another by randomizing the latter. 

Take two probability distributions p = col(pi,p 2 , ■■■,Pd) and q = co\{qi,q 2 , ...,qd), where d can be infinite. We say 
that p majorizes q, and denote it by p q, if and only if 

m m 

Vm<d, (5.1) 

n—1 n—1 

where p'*' and qf are the original vectors with their components rearranged in decreasing order. 

This characterization of the majorization relation is interesting from an operational point of view, since it is easy to 
check numerically if two vectors satisfy this condition. Nevertheless, it can be proven that p q is strictly equivalent 
to two other operational relations: 

(ml) For every concave function h{x), it is satisfied hfpn) < X]n=i 

(m2) q can be obtained from p as q = Dp, where Z? is a column-stochastic matrisQ- 

These relations a very interesting from an interpretational point of view. First, note that the entropy is a concave 
function, and hence, the first relation says that the entropy of q is larger than the entropy of p; now, as we have seen 
more entropy means more information loss, which can be somehow interpreted as more disorder. Second, note that 
any column-stochastic matrix can be written as a convex sum of permutations; hence, the second relation says that 
q can be obtained from p by applying a random mixture of permutations on the latter. Hence, both (ml) and (m2) 
seem to state that q is more disordered than p. 

5.2 Majorization and ensemble decompositions of a state 

As a first simple application of majorization to quantum mechanics, we give an answer to the following question: under 
which conditions can we find an ensemble decomposition based on a probability distribution w, say {wk, \'pk)}k=i, 2 ,...,Ni 
of a density operator p with diagonal representation p = I’’n 

The answer is rather simple: only when A w. Of course, ii N ^ d, zeros are added in the vector with smaller 
dimensionality to match the dimensions of A and w. 

This result is indeed a simple consequence of what we learned in Chapter [5] concerning the freedom to represent a 
given density operator by different ensemble decompositions. In particular, if p can be represented by the ensembles 
{wk, \pk)}k and {A„, |r„)}n, there must exist a left-unitary matrix U for which 

y/ffnWi) = (5.2) 

n 

now, taking the inner product of this expression with itself, and using the orthonormality of the {\rn)}n set, we get 

Wk = J2\^kn\'' Xn- (5.3) 


^ A square matrix is column-stochastic if its elements are real and positive, its columns sum to one, and its rows sum to less than one. 
Most of the literature on the connection between majorization and quantum information studies finite-dimensional systems, in which case 
it can be shown that column-stochastic matrices are also doubly-stochastic (columns and rows both sum to one). One needs the slightly 
more general definition of column-stochastic to cope with infinite dimensional spaces, as we will do in the next chapter. 
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5 . Majorization in quantum mechanics 


Finally, as U is left-unitary, the matrix with elements |Wfcn|^ is column-stochastic. 

5.3 Majorization and the transformation of entangled states 

The next application of majorization theory to quantum mechanics appears when answering another question of 
paramount importance in information theory: given a bipartite state \'>P)ab S ft a ®'Hb, under which conditions can 
it be transformed deterministically into another state \^p)ab G Ha 0 'Hb if Alice and Bob are allowed to use only 
LOGO protocols? 

By ‘deterministically’ we mean that given a complete set of LOGO protocols {fijocc}® transformation \iP)ab — >■ 
\'.p)ab succeeds for all of them, that is, 'S’locc[IV’)ab(V'|] oc |i^)ab(</ 5| Vj. If, for example, the transformation works 
only for fjjocc’ then Alice and Bob will be able to transform |V’)ab into \^p)ab only some of 

the time, in particular with probability po = hence the transformation fails with probability 

(1 — po): that is, Alice and Bob’s strategy works only probabilistically. 

To answer this question, consider the diagonal representations of the reduced states p^ = ti b{I' d)AB{'d\} = 
^tdn)ip{rn\ and p\ = trB{|<p)AB(p|} = Yfn=i^nVri)ip{rn\- Then, Alice and Bob can transform |V’)ab into 
|<p)ab via an LOGG strategy, if and only if in which case we use the symbolic notation \tp)AB -< \'p)ab- 

Note that since the entanglement entropy is a concave function of the eigenvalues of the reduced density operator, 
this majorization relation implies that \'4))ab can only be transformed deterministically via an LOGG protocol into 
states of lower entanglement, that is, E\\d)AB\ > E\\(p)ab]- 

It is also possible to prove that if |V')ab can be transformed into \d)AB deterministically via an LOGG protocol, 
it can always be done with a one-way-direct LOGG protocol of the following form: Bob performs a measurement 
described by some POVM {Hj }and communicates the outcome, say ‘jo’, to Alice, who applies a unitary 
quantum operation Aj^ chosen from a pre-agreed set of unitaries in one to one correspondence with the 

possible outcomes of Bob’s measurement. Hence, if are the measurement operators associated to Bob’s 

POVM, the transformation is accomplished as 

\y})ABO^{Aj®i){i®Bj)\d)AB Vj. (5.4) 



6. QUANTUM INFORMATION WITH CONTINUOUS VARIABLES 


6.1 Infinite dimensional Hilbert spaces: the Harmonic Oscillator 


The harmonic oscillator is the prototypical system which is described quantum mechanically by an infinite dimensional 
Hilbert space. To see this, let’s find the eigenstates of its Hamiltonian, which is given by 



moj 


-9 , 


( 6 . 1 ) 


being m the mass of the oscillator and lu its oscillation frequency. According to the postulates of quantum mechanics, 
the position q and momentum p satisfy the commutation relation 


[q,p]=ih. (6.2) 

We will always work with dimensionless versions of them, the so-called X- and P-quadratures (although we may keep 
using the names ‘position’ and ‘momentum’ most of the time) 

i = “d r = (6-3) 

which satisfy the commutator 

[A,P] = 2i, (6.4) 

and in terms of which the Hamiltonian reads 

R = (6.5) 

In order to find the eigensystem of this operator, we decompose the quadratures as 

X = -I-a, and P = i(d^ — d), (6.6) 

where the operators d and d^, known as the annihilation and creation operators, satisfy the commutation relation 

[d,d'l'] = l; (6.7) 

in terms of these operators, the Hamiltonian is rewritten as 

R = da;(d^d-I-1/2), (6-8) 

and hence the problem has been reduced to finding the eigensystem of the so-called number operator N = cha. 

Let us call n to a generic real number contained in the spectrum of N, whose corresponding eigenvector we denote 
by |n), so that, N\n) = n\n). We normalize the vectors to one by definition, that is, (n|n) = 1 Vn. The eigensystem 
of N is readily found from the following two properties: 

• iV is a positive operator, as for any vector \'ijj) it is satisfied ('i/|iV|^/) = {d\ijj), d\'d)) > 0. When applied to its 
eigenvectors, this property forbids the existence of negative eigenvalues, that is, n > 0. 

• Using the commutation relatior0 [iV,d] = —a, it is trivial to show that the vector d|n) is also an eigenvector of 
N with eigenvalue n — 1. Similarly, from the commutation relation [N, d^^] = d^ it is found that the vector d^|n) 
is an eigenvector of N with eigenvalue n -I- 1. 


^ This is straightforward to find by using the property [AB, C] = A[B, C] + [A, C]B, valid for any three operators A, B, and C. 
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6 . Quantum information with continuous variables 


These two properties imply that the spectrum of N is the set of natural numbers n G {0,1, 2,...} = N, and that the 
eigenvector |0) corresponding to n = 0 must satisfy a|0) = 0; otherwise it would be possible to find negative eigenvalues, 
hence contradicting the positivity of N. Thus, the set of eigenvectors is an infinite, countable set. Moreover, 

using the property a|0) = 0 and the commutation relations, it is easy to prove that the eigenvectors corresponding 
to different eigenvalues are orthogonal, that is, {n\m) = 6nm- Finally, according to the axioms of quantum mechanics 
only the vectors normalized to one are physically relevant. Hence, we conclude that the Hilbert space generated by 
the eigenvectors of N is isomorphic to P (oo)—see Section lA.31 — , and hence it is an infinite-dimensional Hilbert space. 

Summarizing, we have been able to prove that the Hilbert space associated to the one-dimensional harmonic 
oscillator is infinite-dimensional. In the process, we have explicitly built an orthonormal basis of this space by using 
the eigenvectors {|'ra)}„gjsj of the number operator N, with the annihilation and creation operators {a, a'^} allowing us 
to move through this set as 

a\n) = y/n\n — 1)^ and a^\n) = ^/n + l|n -I-1), (6-9) 

the factors in the square roots being easily found from normalization requirements. 

In contrast to the number operator, which has a discrete spectrum, the quadrature operators possess a pure 
continuous spectrum; let us focus on the X operator, whose eigenvectors we denote by {|a:)}x6R with corresponding 
eigenvalues {a:}x6R, that is, 

X\x) = x\x). (6.10) 

In order to prove that X has a pure continuous spectrum, just note that, from the relation 


exp 


yP X exp --yP =X + y, 


( 6 . 11 ) 


which is easily found via the Baker-Campbell-Haussdorf lemmgH, it follows that if |a;) is an eigenvector of X with x 
eigenvalue, then the vector exp(—i2/P/2)|x) is also an eigenvector of X with eigenvalue x + y. Now, as this holds for 
any real y, we conclude that the spectrum of X is the whole real line. Moreover, being a self-adjoint operator, one can 
use the eigensystem X as a continuous basis of the Hilbert space of the oscillator by using the Dirac normalization 
{x\y) = 5{x — y). The same results can be obtained for the P operator, whose eigenvectors we denote by {|p)}pgR 
with corresponding eigenvalues {p}pgR, that is, 

P\p)=p\p). (6.13) 

It is not difficult to prove that there exists a Fourier transform relation between the position and momentum bases, 
that is. 


\P) = 


+ 00 


dx 

j ' 

v47r 


exp -px \x) 


k) = 


^+oo 


dp 


exp -TiPx Ip). 


To this aim we now prove that 


{x\p) = _ exp(ixp/2). 


•\/4^ 


First note that the commutator [X, P] = 2i implies that 

{x\P\x') = 


2i5 {x — x') 


X — X' 


and hence 


(6.14) 

(6.15) 

(6.16) 


(x|P|V') = / ^\ x'\ 


= / dx 


, 2i6 (x — x') 


X — X' 


X — X' 


{x\'tp) + {x' - x) 


d{a 


dx 


E 


{x' — a;)” d”(x|' 
nl dx'^ 


(6.17) 


the order zero of the Taylor expansion is zero because the Kernel is antisymmetric around x, while the terms of order 
two or above give zero as well after integrating them. This means that 


(a;|P|'0) = —2i 


dixl-ip) 

dx ' 


^ This lemma reads 


‘ r) IV V V V 


(6.18) 


( 6 . 12 ) 


and is valid for two general operators A and B. 
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which applied to j'l/') = \p) yields the differential equation 


= -2i^, (6^19) 

which has (16.151) as its solution, the factor l/-\/4^ coming from the Dirac normalization of the |p) vectors. 

As an example of the use of these continuous representations, we now find the position representation of the number 
states, which we write as 

\n) = / dx'ipnix)\x). (6.20) 

Jr 

As a first step we find the projection of vacuum onto a position eigenstate, the so-called ground state wave function 
tpo{x) = (a;|0), from 

0 = (a;|a|0) = ^{x\{X + iP)|0) = + 2^^ i’oix), (6.21) 

where we have used H6.18L which is a differential equation for ifuix) having 

V’o(a;) = ^ 2 ^ y/a exp(-a;V4), (6.22) 

as its solution; the factor (27r)“^/^ is found by imposing the normalization 

(0|0) = / dx-iflix) = 1. (6.23) 


Now, the projection of any number state |n) onto a position eigenstate (the excited wave function) is found from 
the ground state wave function as 


1 


which, reminding the Rodrigues formula for the Hermite polynomials 


fjnix) = {x\n) = ^{x\d^^\0) = 

^/n\ Vn!2’^ 




3 ^- 2 ^) Mx). 


H„{x/V2) = 2 exp(a;^/4) — 2^^ exp(—x^/d), 


(6.24) 


(6.25) 


leads to the simple expression 


V2"+i/27ri/2„! 


(6.26) 


Note that not being vectors contained in the Hilbert space of the oscillator (they cannot be properly normalized), 
the position and momentum eigenvectors cannot correspond to physical states; nevertheless, we will see that they can 
be understood as a (unphysical) limit of some physical states (the squeezed states). 

Finally, let us stress that even though all that we are going to discuss in what follows applies to a general bosonic 
system, that is, a system described by a collection of harmonic oscillators, we will always have in mind the electromag¬ 
netic field (light, in particular), which can be described as a set of modes with well defined polarization, frequency, 
and spatial profile, each of which behaves as the mechanical harmonic oscillator that we have introduced. 


6.2 The quantum harmonic oscillator in phase space I: The Wigner function 

As the position and momentum do not have common eigenstates, and moreover, their eigenstates cannot correspond to 
physical states of the oscillator, one concludes that these observables cannot take definite values in quantum mechanics; 
given the state p, the best one can offer is the probability density function which will dictate the statistics of a 
measurement of these observables, that is, P{x) = (a;|p|a:) and P{p) = {p\p\p)- In other words, quantum mechanically 
there are not well defined trajectories in phase space, the position and momentum of the oscillator are always affected 
by some {quantum) noise. 

It follows naturally the following question; is it then possible to describe quantum mechanics as a probability 
distribution defined in phase space which blurs the classical trajectories? As we are about to see, the answer is only 
partially positive, as quantum noise is a lot much subtle than common classical noise. 
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6 . Quantum information with continuous variables 


A logical way of building such a phase space distribution, say Wp{x,p), is as the one having the probability density 
functions P{x) and P(p) as its marginals, that is, 

P{x) = f dpWp{x,p), and P{p) = f dxWp{x,p). (6.27) 

Jk Jk. 

It is possible to show that this distribution is uniquely defined by [S] 

Wp{x,p) = ^ j dyexp {x + y/‘2\p\x - y/2), (6.28) 

which is known as the Wigner function. It is immediate to check that this distribution has the proper marginals, and 
that it is normalized, i.e., dxdpWp{x,p) = 1. The proof of its uniqueness is not that simple though. 

As we will see in a moment, given the state p of the oscillator, the quantum expectation value of an observable 
0{X,P) can be evaluated via the Wigner function as 

{d{X,P)) =ti{pd} = [ dxdpW^{x,p)0^^\x,p), (6.29) 

dR2 


where {x,p) is obtained by writing d(X, P) as a symmetric function of X and P with the help of the commutation 
relation (16.41) , and changing these operators by the real variables x and p, respectively. 

Taking into account that the prescription to find the quantum operator associated to a classical observable 0{x,p) 
consists precisely in symmetrize it with respect to x and p, and then change the position and momentum by the corre¬ 
sponding self-adjoint operators (what guarantees the self-adjointness of the remaining operator), this result reinforces 
the interpretation of Wp{x,p) as a probability distribution in phase space, and hence of quantum mechanics as noise 
acting onto the classical trajectories. However, we will find out later that this distribution can be negative, and hence 
it is not a true 2D probability density function. Note that, actually, if Wp(a;,p) was positive all over phase space for 
any p, quantum mechanics could be simulated with classical noise, while quantum mechanics has been proven itself 
to go well beyond classical mechanics in many experiments. 

Given a state p of the oscillator, expression (I6.28|) allows us to compute the corresponding Wigner function. 
However, there is a much more convenient way of writing the Wigner distribution, which is based in the so-called 
displacement operator 


i)(r) = exp 


-R^ Hr 


(6.30) 


where r = col(a;,p) is the coordinate vector in phase space (which is usually denoted by the displacement of the 

is known as the symplectic form. 

Note that with this matrix notation the position-momentum commutators read [Rj,Ri] = 2iQji. We will learn 
more physical things about the displacement operator in the following sections, but for now, just take it as a useful 
mathematical object. 

Note that using the formula 


operator), R = col(A,P) is the corresponding vector operator, and H = 


exp(A + B) = exp(—[A, B]/2) exp(A) exp(.B), 


(6.31) 


valid for operators A and B which commute with their commutator, this operator can be written as a concatenation 
of two individual momentum and position translations plus some phase: 



i 


i 


i A 

exp 

1 

exp 

^pX 

exp 

- tP 

[ 2^^ \ 


(6.32) 


which allows us to write 


tr{D{r)} = / dx'dp'{x'\e^P^/^e-'^^^/^\p'){p'\x') 

<^Jpx'/2 [ (V iV/2 
4^® ./„4^® ’ 




S{p) 


S{x) 


(6.33) 

(6.34) 


arriving to the identity 


tr{Z)(r)} = 47r5^^^(r), 


( 6 . 35 ) 
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which will be useful in many situations. Another important property, trivially proved from (I6.31L is 


D{r)D{r') = exp 


—-r^idr' ) D{r + r'), 


(6.36) 


and therefore, except for a phase, the concatenation of two displacement operators is equivalent to a single displacement 
operator with the sum of the displacements. Note that the phase is zero only when the condition xp' = px' is met, 
although it plays no physical role when applied to a state of the system. 

The first step in order to find the Wigner function of a given state p is to define the characteristic function 

Xp{r) =ti{D{r)p} p= f ^Xp(r)i)1'(r), (6.37) 

jR2 dTT 

and then the Wigner function is obtained as its Fourier transform 

^ Xp(s) = / d^rWp(r)e-^^"^- (6.38) 

dR2 (47r) Jr2 

it is not difficult to show that this alternative definition of the Wigner function leads to the original one given by 
(I6.28|l . However, this definition simplifies a lot many derivations. For example, from (16.3711 and (I6.38L it is immediate 
to prove that 

trip} = 1 ^ Xp(0) = 1 ^ [ (PrWp{r) = 1, (6.39) 

that is, the normalization of the Wigner function. Moreover, evaluating the trace of ()6.37l) in the position eigenbasis, 
we get 

Xp(0,p) = J dx{x\e^^^p\x) = j dxe^'^^/‘^{x\p\x), (6.40) 

which using the right hand side of (I6.38P directly implies that 


{x\p\x) 



(6.41) 


that is, the Wigner function has the position probability density function as one of its marginals. The other marginal 
is the momentum probability density function, as is proved from Xp(a;,0) in a similar fashion. 

Going one step further, using this definitions it is actually quite simple to prove (16.291) . In particular, it is easy to 
check that the expectation value of the symmetrically ordered product can be obtain from the characteristic 


function as 

P)j+i 

{(X>/>•)'■>) = (-)y2i)«‘ g^X,{r) 

(6.42) 

r^O 

leading to 

((A^A)(")) = [ dxdpWp{r)x^p\ 

(6.43) 


after using (I6.37L in agreement with (I6.29|) . 

In general we will not deal with a single harmonic oscillator (a single mode of the light field), but with a collection 
of, say, N harmonic oscillators {N modes of light). Let us define again the coordinate vector in the complete phase 
space by r = (xi,pi,X 2 ,P 2 , ■■■,xntPn), and the corresponding vector operator 

R = col(Ai, A, ^ 2 , A, Xn, Pn). (6.44) 


in terms of which the commutation relations can be rewritten as 


[Rj,Ri] = 2i{HN)ji, 


(6.45) 


where 


N 


— 


n 


n 


n 


(= 


(6.46) 
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6 . Quantum information with continuous variables 


is the symplectic form of N modes (in what follows we will suppress the subindex N except when needed). The state 
of the system acts now onto the tensor product of the Hilbert spaces of the modes, and so does the displacement 
operator, which is now defined as 


D{r) = exp 


-R^ Hr 


= D{ri) (g) D{r 2 ) g)... g T*(rAr), 


and satisfied tr{£>(r)} = as well as (16.361) : the characteristic function is defined as before 


Xp{r) = tr{D{r)p} 


P = 


RSN- 


(47r) 


N 


Xp(r)£)1'(r), 


and the Wigner function as its 2iV-dimensional Fourier transform 

(P^S 


Wp(r) = 


/R2W 


(47r) 


2N 


Xp(s)e 


r2r 


X^(s) = / P^rWpir)e-^ 


— Qr 


(6.47) 


(6.49) 


(6.50) 


We saw in the previous chapters that there is an operation that plays an important role when dealing with composite 
Hilbert spaces: the partial trace. For example, in the case of the N harmonic oscillators being in a state p, imagine 
that we want to trace out the last one, obtaining the reduced state of the remaining — 1 oscillators, pn = trjv{;o}; 
let’s see what this means in phase space. From (16.371) we see that the characteristic function of the reduced state pn 
is just the original one with the phase space coordinates of the traced oscillator set to zero, that is, 


XR(r{jv-i}) = tr{b{r{N-i})PR} = Xp(i’{Ar-i}, 0), 


(6.51) 


where we use the notation = col(ri, r 2 ,..., r 7 v_i), and therefore, the Wigner function associated to the reduced 

state pn can be found by integrating out the phase-space variables of the corresponding oscillator, that is, 

WR(r("-i)) = = J d2r^Wp(r{^-i},r^). (6.52) 

Let us end this section with an interesting example of a Wigner function: that of a number state p = |n)(n|. 
Although not straightforwardly, it possible to show that the corresponding Wigner function is given by 


W\n){,x,p) = -f p2)exp 

where Ln{z) is the Laguerre polynomial of order n, which can be found from the Rodrigues formula 

exp( 2 :) P 


Lnib — 


i! dz' 


[z"exp(-z)]. 


(6.53) 


(6.54) 


Note that, as commented above, this function has negative regions (except for n = 0, corresponding to the vacuum 
state), and therefore, it cannot be simulated with any source of classical noise. For example, for odd n it is always 
negative at the origin of phase space, since Ln(0) = 1 Vn. 


6.3 Gaussian continuous variables systems 

6.3.1 Gaussian states 
General definition 

A particularly important class of quantum mechanical states of the harmonic oscillator are the so-called Gaussian 
states., that is, states which have a Gaussian Wigner function; as we will see, these are the type of light states which 
are naturally generated in the laboratory, although we will also show how to design experimental schemes whose 
purpose is the generation of non-Gaussian states. 


® Note tliat the 2Af-dimensional Dirac delta function can be written as 


Jm 


I{2W (47r)2^ 


(6.48) 
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The Wigner function of an 


arbitrary Gaussian state has then the forrrQ 


W{r) 


1 

27r-\/det V 


exp 






(6.56) 


where we have defined the mean vector 

r = (R)=col ((1),(P)), 

and the covariance matrix 

r (51^) U{^x,6P}) 

msx,sp}) {SP^) . 


whose elements are given by 


(6.57) 

(6.58) 


(6.59) 


in this expression we have used the notation 5A = A — {A), and denoted the anticommutator by curly brackets. Note 
that number states do not belong to these class of states (save vacuum, as we show below), but we will show that 
many other several interesting states do. Note also that Gaussian states are completely defined by their first and 
second moments, which is why we will sometimes denote by po{f,V) a given Gaussian state. Finally, we would like 
to remark that the mean photon number of the state is given by 


W = + if» - i 


(6.60) 


Classically, any covariance matrix is allowed, as long as it is real, symmetric, and positive definite. A quantum 
mechanical harmonic oscillator has the added constrain det F > 1, what comes from the uncertainty principle of 
position and momentum; indeed, the following proof is quite reminiscent of the proof of the uncertainty principle: 

detv = {5X^){5P^)-\{{5X,5P}p > \{5X5P)\'^-\{{5X,5P}p (6.61) 

4 Schwartz 4 

= 1\{[SX,SP]) + {{SX,SP})\^ - \{{SX,SP}p = 1. (6.62) 

imaginary real 


A Gaussian state corresponding to any real, symmetric covariance matrix satisfying this condition is physically achiev¬ 
able. Note that this condition directly implies the positivity of the covariance matrix. 

In the case of dealing with a Gaussian state of N modes, their Wigner function takes the following form in the 
complete phase space: 


W{r) 


1 

- , exp 

(27r)^VditF 






(6.63) 


where now r is a vector with 2N components Vj = (Rj), and F is a 2N x 2N matrix with elements Vji = {{SRjSRi}) /2. 
Note that the mean photon number is given in this multi-mode case by 


Ew) = ■ y- 


(6.64) 


In this case, the necessary and sufficient condition for a real, symmetric 2N x 2N matrix F to correspond to a 
physical quantum state of N oscillators is that [5] 


F-biG>0, (6.65) 

where the inequality must be understood as “positive semidefinite”, what is again linked to the uncertainty principle, 
and implies the positivity of F. It is possible to prove that any Gaussian state for which the eigenvalues of F -|- ifl 

^ When dealing with Gaussian states, the following integral is quite useful: 

d^^rexp ^ exp , 

where r and A is a nonsingular 2N X 2N matrix. 



(6.55) 
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are zero (equivalently, which saturates the uncertainty principle) is pure, and there exists no other non-Gaussian state 
with the same covariance matrix. 

In many situations it is interesting to understand the state of the N oscillators as a bipartite state of M oscillators 
plus another M' oscillators (M + M' = N), in which case we talk about an M x M' continuous variable system. 
Consider a Gaussian state /5G(r, V) of the N oscillators, whose mean vector and covariance matrix we write as 


r = col(x,x'). 


and 



C 

W' ’ 


( 6 . 66 ) 


where x € x' € , W and W real, symmetric matrices of dimensions 2M x 2M and 2M' x 2M', respectively, 

and C a real matrix of dimensions 2M x 2M'; then, it is easy to prove that the state of the first M modes after tracing 
out the remaining M' modes, is the Gaussian state 


/5g(x, W) = trM+i,M+2,...,M+M'{/3G(r, G)}, 


(6.67) 


that is, tracing out a mode in a Gaussian state is equivalent to remove its corresponding entries in the mean vector, 
as well as its rows and columns in the covariance matrix. In order to prove this we make use of the characteristic 
function, which for the general Gaussian Wigner function (I6.6dl) takes the form 


X(s) = exp 


+ —r^Vls 
8 4 


( 6 . 68 ) 


then, by substituting (16.661) in this expression, and remembering that tracing out any mode is equivalent to setting to 
zero the corresponding phase space variables in the characteristic function, we can write the characteristic function of 
the reduced state of the first M modes as 


Xfl(s{M}) = exp 




(6.69) 


what proves (16.671) . 

Note hnally that when the states of the partitions are uncorrelated, that is, 

PG(r, V) = /3g(x, W') 0/5G(x',kk')) 

the covariance matrix of the Wigner function (16.631) can be written as a direct sum 

V = W(BW'. 


(6.70) 


(6.71) 


Examples of Gaussian states 

The vacuum state. As commented above, number states are not Gaussian; nevertheless, there is one exception: the 
vacuum state |0). To see this, just note that its Wigner function (16.531) can be written as 


W\o){x,p) = ^exp ( - 

that is, a Gaussian distribution like (16.561) with 

r = 0, and V = 




I 0 
0 1 


= X, 


■2x2- 


(6.72) 


(6.73) 


Hence, in Gaussian notation |0)(0| = /5 g(0,I2x2)- 

Thermal states. As explained in Chapter [21 the mixedness of the state of a system can be related to the amount 
of information which has been lost to another system inaccessible to us, that is, to the amount of correlations shared 
with this second system. Given a system whose associated Hilbert space has dimension d, and is spanned by some 
orthonormal basis {|j)}j=i, 2 ,...,d, we already saw that its maximally mixed state is 
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where Id is the identity of the c?-dimensional Hilbert space. Being proportional to the identity, this state is invariant 
under changes of basis, and hence, the eigenvalues of any observable of the system are equally likely; this is in 
concordance with what one expects intuitively from a state which has leaked the maximum amount of information to 
another system. 

For infinite dimensional Hilbert spaces {d —>■ oo) this state is not physical since it has infinite energy, that is, 
tr{/3fV} —>■ oo. Hence, the question of which is the maximally mixed state in infinite dimension makes sense only if 
one adds the energy constraint tr{plV} = h, where h is a positive real. It is possible to show that the state which 
maximizes the von Neumann entropy subject to this constraint is 

oo 

/5th(h) = ^ ^ i+„ |n)(?t|, (6.75) 

n=o (1 + 

which is still diagonal in the number state basis, but has not a flat distribution for the number of photons. Interestingly, 
it can be appreciated that this distribution is the one expected for a collection of bosons at thermal equilibrium with 
an environment, the Bose-Einstein distribution. This is why this state is known as the thermal state, whose von 
Neumann entropy reads 

‘S'[/5th(h)] = (h + 1) log(fi + 1) - hlogh = Sthih). (6.76) 

It is not difficult to see that this state is Gaussian (later on we will actually prove it by simple means), and that 
it is defined by a zero mean vector, and a covariance matrix 

Hth(h) = (2h +1)12x2, (6.77) 

that is, pth(ji) = /5g[ 0, Hh(h)]. Note that the vacuum state can be seen as a thermal state with zero mean photon 
number. 

In the next section we will learn that any A^-mode Gaussian state can be decomposed into N uncorrelated thermal 
states (Williamson’s theorem), and hence thermal states can be seen as the most fundamental Gaussian states. 

6.3.2 Gaussian unitaries 
General deSnition 

Gonsider a unitary transformation U = exp{H/ih) acting on the state of the light field. We say that such a unitary is 
Gaussian, if it maps Gaussian states into Gaussian states. 

Let us define the vector operators 


a = col(ai, 02 ,..., Oat), and a"*" = col(a|, o^,..., a]^); (6.78) 

it is quite intuitive that any Gaussian unitary will come from a Hamiltonian having only linear or bilinear terms, that 

is, 

H = ih(^a. ■ + a'^'^C/a’^) + H.c., (6.79) 

for some vector a £ C'^, and some symmetric, complex N x N matrices E and Q. Let’s try to understand why 
this is so by analyzing the physical meaning of each term. The first term corresponds to the injection of photons 
in the modes; the second term comprises all the energy shifts Ojdj, as well as the creation of a photon in one mode 
via the annihilation of a photon of another mode, Ojdi; the last term takes into account the possibility of generating 
two photons simultaneously, In other words, all the possible one-body and two-body interactions are taken into 

account in this Hamiltonian. It is obvious that if we want to transform Gaussian states into Gaussian states, the 
Hamiltonian cannot contain many-body interactions beyond these, because otherwise one would create correlations 
which go beyond the ones captured by the covariance matrix, and hence the final state would not be completely 
characterized by its first and second moments. 

In the Heisenberg picture, a Gaussian unitary would then induce a so-called Bogoliubov transformation 

a—>■ E^alJ = Aa +Ba^ + a, (6.80) 

where the form of the complex N x N matrices A and B in terms of (a,iF,Q) is unimportant for our purposes; the 
only restriction on these is that they have to satisfy AB'’" = BA^ and AA"^ = BB^ -|-lArxAr, in order to preserve the 
commutation relations of the creation and annihilation operators (Inxn is the identity matrix of dimension N). 
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Instead of writing this linear transformation in terms of the boson operators, one can write it in terms of the 
position and momenta, or more compactly, in terms of the R vector operator: 


R ^ U^RU = 5R + d, 


(6.81) 


where, once again, the dependence of d G and the real 2N x 2N matrix S in the previous transformation 

parameters is unimportant for our purposes; the only relevant thing is that, in order to preserve the commutation 
relations of the quadratures, S has to satisfy 

sns^ = n, (6.82) 

that is, S must be a symplectic matrix. In the following we will denote Gaussian unitaries by UG{d,S) to stress the 
fact that they are completely characterized by d and S. 

The transformation induced onto the system by the Gaussian unitary Uc^d, S) is easily described in the Schrddinger 
picture as well if the states are represented by the Wigner function. To see this, let us first note that the displacement 
operator is transformed by the action of this unitary as 


17^^(r)17G = exp 

^(5R-kd)^fIr 

exp 

^{R^ns-^ + d^n)r 

= b{S ^r)exp 

-d^flr 

2 








(6.83) 


where the identity = QS ^ follows from (16.821) and (16.461) . Therefore, given the initial state p of the system, the 
corresponding characteristic function is transformed as 

XOopf/t (r) = tr{b{v)UGpUl} = e-^'^"^^ir{b{S-^v)p} = (6.84) 

and the Wigner function as 


W, 


UgpUI 


(r) = 


/ 

jR2Ar 




(47r) 


2N 






Xp{S s)e 




/R2Ar 


(47r) 


2N 


Xp(z)f 


d) 


^ Wp[S~br-d)]. 


(6.85) 

In the case of Gaussian states the situation is even more simple: one only needs to hnd how the transformation 
affects the first and second moments of the state. Using (I6.81L it is straightforward to show that the transformation 
induced by the Gaussian unitary UG{d,S) on the mean vector r and the covariance matrix V of any state is 


5r + d. 


and V 


SVS^ 


( 6 . 86 ) 


Taking again the M x M' partition of the multi-mode system, note that when the unitary transformation acts 
independently on each partition, that is, 

UGid,S) = Ug{1,C) ® Ug{1',C'), (6.87) 

where 1 G 1' G , and C and £' are symplectic matrices of dimensions 2M x 2M and 2M' x 2M', respectively, 
we can write 

d =(!,!'), and 5 = £©£'. (6.88) 

Finally, we would like to remark that a Gaussian unitary transformation is called passive when it conserves the 
number of mean photons active if it changes it. Now, given the transformation (I6.86|) . and reminding 

that the mean number of total photons is proportional to the square modulus of the mean vector, r^, and the trace of 
the covariance matrix, trU, see (16.6411 . a Gaussian unitary will be passive if and only if 

d = 0, and S^S = I 2 Nx 2 N, (6.89) 

the second condition meaning that its associated symplectic transformation must be orthogonal, that is, 5^ = 


Examples of Gaussian unitaries and more Gaussian states 
The displacement operator and coherent states. Consider the unitary operator 

b{a) = exp(aa^ — a*a), (6.90) 

which, using the formula (16.311) can be written in the following two equivalent ways 

(a) = exp (—|ap/2) exp (aa’^) exp (—a*a), or (a) = exp (|q;P/ 2) exp (—a*a) exp (aa^) , (6.91) 
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which we will refer to as its normal and antinormal forms, respectively. 

Using the Baker-Campbell-Haussdorf lemma (I6.12|l . it is fairly simple to prove that this operator transforms the 
annihilation operator as 

a ^ {a)aD{a) = a + a, (6.92) 

or, in terms of the quadratures 

X ^ {a)XD{a) = X + Xa, and P ^ D\a)PD{a) = P + Pa, (6.93) 

where a = {xa + \pa)l"2 with Xq, G K and Pa G R. This unitary operator is then called the displacement operator 
because it allows us to perform translations in phase space; indeed, it is exactly the operator that we defined in the 
previous section, see (I6.30L what is easily shown by rewriting (I6.90|) in terms of the position and momentum operators. 
As a Gaussian unitary we have D[a) = UG{da,l2x2) with 

d„ = col(a;a,Pa)- (6.94) 

The states obtained by displacing the vacuum state are known as coherent states. Using the normal form 
of the displacement operator, it is easy to obtain 


|a) 


n I 


n—0 


(6.95) 


They are Gaussian states with the same covariance matrix as vacuum, but with a non-zero mean vector, that is. 


r = d„, and U=J 2 x 2 , 


(6.96) 


or in Gaussian notation, |a)(Q!| = pc (da, 12 x 2 )- Hence, these states have the same noise properties as vacuum, but 
describe a bright light beam, that is, a light beam with non-zero mean field. In fact, they are a fair approximation to 
the state describing the beam coming out from a (phase-locked) laser. 

From a mathematical point of view, they are the eigenstates of the annihilation operator, that is, a\a) = a\a). 
Later on we will learn that, even though the annihilation operator is not self-adjoint, we can build a POVM-based 
measurement which has its eigenvalues as the possible outcomes [heterodyne detection). 

The squeezing operator and squeezed states. Gonsider now the squeezing operator 

S{r) = expi^a^ - (6.97) 


where r G [0,oo[. This operator is implemented experimentally for an optical mode of frequency ojq by pumping with 
a strong laser beam of twice that frequency a crystal with second order nonlinearity; pairs photons of frequency ojq 
are generated via the so-called spontaneous parametric down-conversion process. 

Using the Baker-Campbell-Haussdorf lemma (16.1211 . it is again simple to prove that this operator transforms the 
annihilation operator as 

a —>■ S\r)dS[r) = a cosh r — a’^sinhr, (6.98) 

or, in terms of the quadratures 

X ^ S\r)XS{r) = e-''X, and P ^ S'f [r)PS[r) = e''P, (6.99) 

so that it is characterized as a Gaussian unitary by S[r) = C/g[0, Q[r)] with 


Q[r) = 




( 6 . 100 ) 


Applying the squeezing operator to a vacuum state, one obtains a so-called squeezed vacuum state. In the number 
state basis, this state is characterized by containing only an even number states, what comes from the fact that the 
squeezing operator generates photons in pairs; its explicit representation in this basis is found to be m 


\r) = 5(010) = 


^ 2"n! 

n—0 


(2n)! 

coshr 


tanh" r|2n). 


( 6 . 101 ) 
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This Gaussian state has zero mean, and covariance matrix 

Vsq{r) = Q(r)Q^(r) = 


( 6 . 102 ) 


that is, |r)(r| = /5g[0, VijqC?’)]- Note that in the limit r —>■ oo the variance of the position goes to zero, while the variance 
of the momentum goes to infinity, and hence in the limit of infinite squeezing the state (jB.lOlIl is an eigenstate of the 
position operator. Note, however, that this limit is unphysical, as the number of photons (N) = sinh^ r diverges, and 
hence, an infinite amount of energy is needed to generate a position eigenstate. 

The phase shift operator. The free evolution of an oscillator (corresponding to the free propagation of an 
optical mode through a linear medium) induces the unitary transformation 


R{9) = exp(—i^a^a). 


known as the phase shift operator, which transforms the annihilation operator as 

a —>■ R\0)aR{O) = exp(i0)a. 


(6.103) 


(6.104) 


or in terms of the position and momentum 

X ^ R^ {6)XR{9) = X cosd + Psm9, and P ^ Rp9)PR{9) = P cos9 — X sm9. 
Hence, as a Gaussian unitary this transformation is characterized by P{9) = Ug[Q,TZ{9)], where 


n{9) 


cos 9 sin 9 
— sin 9 cos 9 ’ 


(6.105) 


(6.106) 


what shows that a phase shift is equivalent to a proper rotation in phase space. 

Note that number states are invariant under this transformation—they are eigenstates of R{9) —, and hence, 
thermal states are invariant under rotations in phase space. This is not the case for coherent or squeezed states. 

The t-wo-mode squeezing operator and t-wo-mode squeezed states. All the unitaries considered so far act 
on a single mode of the electromagnetic field, and hence, they cannot be used to induce entanglement between several 
modes. In this example we consider the two-mode squeezing operator 

Si 2 {r) = exp(raia 2 — ra\a\), (6.107) 

which can be implemented experimentally via a nonlinear crystal as the squeezing operator (16.971) . but now in a regime 
in which the down-converted photons are distinguishable either in frequency, and/or polarization, and/or spatial mode. 
Under the action of this operator, the annihilation operators are transformed as 


fii —>■ .Sj 2 (''’)®i‘^i 2 (^) = Si coshr — aj sinhr, (6.108) 

0.2 —t Sl 2 {f)a 2 Si 2 {r) = 02 coshr — oj sinhr, (6.109) 

or in terms of the quadratures 

Xi —>• 5j2(?')-Ai5i2(r) = Xi coshr — X 2 sinhr, (6.110) 

Pi Sl2{r)PiSi2{r) = A coshr-I-P 2 sinhr, (6.111) 

X 2 —t A2(^)^2<5i2(?') = X 2 coshr — Xi sinhr, (6.112) 

P 2 ^ Sl 2 {r)P 2 Si 2 {r) = A coshr-I-A sinhr. (6.113) 


Hence, as a Gaussian unitary this transformation is characterized by 5'i2(r) 


Gg[ 0, Qi 2 (r)] with 


212(1-) 


2 ^ 2 x 2 coshr —2^sinhr 
—2^sinhr X 2 x 2 COshr ’ 


(6.114) 


where Z = diag(l, —1). 

Applying the two-mode squeezing operator to a vacuum state, one obtains a so-called two-mode squeezed vacuum 
state. In the number state basis, this state is characterized by a perfectly correlated statistics of the number of quanta 
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in the modes, what again comes from the fact that the two-mode squeezing operator generates photons in pairs; its 
explicit representation in this basis is m 


1 °° 

'-}i2 = ^i2(?')|0, 0) =- - — y^tanh^rln,?^) 

coshr 


(6.115) 


n=0 


where we have used the notation \n) ® \m) = |n,TO). This Gaussian state has zero mean, and covariance matrix 


V'2sq(r) = Qi2{r)Ql2{r) 


l 2 x 2 C 0 sh 2 ?’ — Zsinh2r 

—Zsinh2r l 2 x 2 COsh 2 r ’ 


(6.116) 


that is, |r)i 2 (r| = ;o[0, V' 2 sq(r-)]. 

Note that by taking the partial trace onto any of its two modes, the two-mode squeezed vacuum state becomes a 
thermal state with mean photon number n = sinh^ r, that is tr 2 {|r)i 2 (r|} = tri{|r)i 2 (r|} = /5th(sinh^r), and hence 
the two-mode squeezed vacuum state can be seen as the purification of a thermal state, what in addition shows that 
it is the maximally entangled state in infinite dimension for a fixed energy. We will come back to the entanglement 
properties of the two-mode squeezed vacuum state in Section [6.3.41 

The beam splitter operator. We are going to analyze only one more type of two-mode unitary transformations, 
the one induced by the so-called beam splitter operator 


Bi 2 (/3) = exp (/3dia2 - I3a\a-^ , (6.117a) 

which can be implemented experimentally by, for example, mixing two optical beams in a beam splitter of transmissivity 
T = cos^ j3. 

Under the action of this operator, the annihilation operators are transformed as 


oi -)■ {P) a^B{p) = oi cos/3-I -02 sin/3 (6.118) 

02 ^ (/3) a2B (/3) = 02 cos/3 — di sin/3. (6.119) 

or in terms of the quadratures 

^ B\^ (/3) XiBi 2 (/3) = cos/3 - X 2 sin/3, (6.120) 

A ^ A 2 (/3) hBi 2 (/3) = A cosP - A sin/3, (6.121) 

X 2 ^ (/3) ^ 2^12 (/3) = X 2 cos/3 -h sin/3, (6.122) 

A ^ A 2 (A ASi 2 (/3) = A cos/3 -f A sin/3. (6.123) 


Hence, as a Gaussian unitary this transformation is characterized by .Bi2(/3) = 17 g[0iS i2(/3)] with 


^12 (A 


12x2 cos/3 -l2x2sin/3 
12x2 sin/3 I 2 x 2 cosh/3 


(6.124) 


It is interesting to note when the states of both modes are coherent, they keep being coherent after the action of 
the beam splitter transformation, as Hi2(/3)H|(^(/3) = 14 x 4 - As an example consider the state |a) 0 |0) (one mode in 
an arbitrary coherent state, and the other in vacuum), which has the Gaussian representation /3G(d,l4x4) with 

d = col(2 Re{a}, 2 lm{a}, 0,0); (6.125) 

after the action of the beam splitter operator, it becomes /5Q(d',l4x4) with 

d' = col(2 Re{a} cos/3, 2 lm{a} cos/3, 2 Re{a} sin/3, 2 lm{a} sin/3), (6.126) 

which is the tensor product of two coherent states, in particular, |acos/3) 0 |asin,S). This is exactly what one expects 
when a laser held is sent through a beam splitter: part of the laser is transmitted, and part is reflected. 
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6.3.3 General Gaussian unitaries and states 

In this section we will use symplectic analysis (or better, ‘symplectic tricks’), to find interesting facts about general 
Gaussian unitary transformations and Gaussian states. 

It is well known in symplectic analysis that any 27V x 27V symplectic matrix S can be decomposed as 


S = /C 


02(g) 


7=1 




(6.127) 


where K, and C are orthogonal, symplectic matrices (this is known as the Euler decomposition of a symplectic matrix, 
or as its Bloch-Messiah reduction). Physically, this means that a general TV-mode unitary transformation can be seen 
as the concatenation of three operations: an 7V-port interferometer mixing all the modetH, TV single-mode squeezers 
acting independently on each mode, and a second 7V-port interferometer. 

As an important example involving two modes, note that the two-mode squeezing transformation can be written 
as 

Qi 2 (r) = Si 2 (-J) [Q(-r) © Q(r)] (6.128) 

what in the Hilbert space means that two-mode squeezed vacuum state can be obtained by mixing a position squeezed 
state with a momentum squeezed state in a 50/50 beam splitter, that is, 

I?')i 2 = 7?i2(-7r/4)[| - r) (g) |r)]; (6.129) 

note that the first beam splitter disappears because the two-mode vacuum state |0) © |0) is invariant under passive 
transformations. Squeezed beams being easily generated in the laboratory, the relation (I6.129|) has allowed us to 
achieve two-mode squeezed beams with high degree of entanglement. 

As a second example, note that for one mode the only passive transformations are the rotations in phase space, 
what means that an arbitrary single-mode Gaussian unitary can be written as the concatenation of a phase shift, a 
squeezing operation, a second phase shift, and a final displacement, that is, 

C/G(a, 6», r, (/) = D{a)R{e)S{r)R{(j)). (6.130) 


Now, it is quite intuitive (although it is not trivial to prove) that any Gaussian state pc having von Neumann entropy 
So can be obtained by applying a unitary transformation onto the thermal state pth(ho) with that same entropy— 
5'th(ho) = So —, that is, 

PGia,0,r,no) = D{a)R{0)S{r)pthino)S^ ir)R\0)D\a), (6.131) 

what means that the covariance matrix of any single-mode Gaussian state can always be decomposed as 


V{0,r,no) = (2ho + l)7^(6l)Q(2r)7^^(6>) = (2ho + 1) 


cosh 2r — cos 20 sinh 2r sin 20 sinh 2r 

sin 20 sinh 2r cosh 2r + cos 20 sinh 2r ’ 


(6.132) 


note that the first phase shift has disappeared because thermal states are invariant under such transformations. 

A second interesting theorem is that of Williamson's, which states that any positive 27V x 27V symmetric matrix 
V can be brought to its diagonal form V® by a symplectic transformation W, that is, 

N 

H = with V®=^vjl 2 ^ 2 - (6.133) 

7=1 

This theorem has a huge application in the world of Gaussian states. Note that, physically, V® can be seen as the 
covariance matrix of TV independent modes in a thermal state with mean photon numbers {hj = (vj — l)/2}j=i^2,...,Ar, 
while the symplectic transformation W corresponds to a Gaussian unitary transformation. Williamson’s theorem is 
then completely equivalent to state that any TV-mode Gaussian state pg(j, V) can be obtained as 


r N 

PG{f, H) = 77G(r, W) j 0/5th[(G " l)/2] 




(6.134) 


^ In optics, an interferometer is just a collection of beam splitters which mix optical beams entering through its input ports. They 
correspond to the most general pasive Gaussian unitary, and are described by a concatenation of single-mode phase shifts and two-mode 
beam splitters. 
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The set {vj}j=i, 2 ,...,N is called the symplectic spectrum of V, so that each i/j is a symplectic eigenvalue, it is possible 
to show that the symplectic spectrum of V can be computed as the absolute values of the eigenvalues of the Hermitian 
matrix iflt/. 

This decomposition is very important, since it allows us to write many properties of Gaussian states and covariance 
matrices in an easy form. For example, the condition (16.6511 which ensures that V is the covariance matrix of a physical 
Gaussian state can be rewritten as 

V >0, and nj > 1 Vj. (6.135) 

As a second important example, note that as unitary transformations do not change the von Neumann entropy, the 
entropy of V) can be directly computed as the sum of the entropies of the corresponding thermal states, that is, 


N N 

5'[pG(f,f^)] = - l)/2] = 

j=i i=i 


where we have defined the function 


aix) = 


log 


X — 1 


log 


a; — 1 


(6.136) 


(6.137) 


which is positive and monotonically increasing for a; > 1. 

It is particularly simple to evaluate the symplectic eigenvalues in the case of dealing with one or two modes. In the 
case of one mode, the trick is to realize that the determinant of the covariance matrix V is invariant under symplectic 
transformations, and hence the sole symplectic eigenvalue reads in this case 

V = V det V. 


For two modes, let us write the covariance matrix in the block form 

V = 


A C 
B 


(6.138) 

(6.139) 


where A = , B = B"^, and C are 2x2 real matrices. In this case, there is an extra symplectic invariant, namely 

A(y) = det A + deti? + 2detC', and hence the symplectic eigenvalues of a general two-mode Gaussian state can be 
obtained from 

det V = and A(V) = 1^1 + 12 ^, (6.140) 

leading to 


^2 _ A(y)± v^A2(F)-4dety 

In terms of the two-mode symplectic invariants, the second condition in (I6.135P is rewritten as 

det V > 1, and A < 1 -|- det V. 


(6.141) 


(6.142) 


It is particularly relevant the case in which the covariance matrix of the two-mode Gaussian state is in the so-called 
standard form 


V = 


a 

0 

Cl 

O' 

0 

a 

0 

C2 

Cl 

0 

h 

0 

0 

C2 

0 

h 


(6.143) 


indeed, it is possible to show that the covariance matrix of any bipartite Gaussian state can be brought to this standard 
form via a local Gaussian unitary transformation Gq = Ug(0, Si)0(7g(O, S 2 ). In the particular case Ci = —C 2 = c > 0, 
the symplectic eigenvalues read 


iy± = 


\/(a + b)^ — 4c^ ± (b — a) 


and the symplectic matrix W satisfying V = can be explicitly found as 

W = 

with 


u;+X2x2 u}—Z 
10—Z W_|_X2x2 


2 1 


1/(0 -I- 6)2 — 4c2 


± 1 


(6.144) 


(6.145) 


(6.146) 
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6.3.4 Gaussian bipartite states and Gaussian entanglement 

In Chapter [3] we introduced the concept of entanglement as correlations between two systems A and B which go 
beyond the ones allowed classically. In this section we particularize those ideas to the case of Gaussian continuous 
variable states. In the following we consider only two modes, that is, a 1 x 1 continuous variable system, although we 
will talk at the end a little about general N x M systems. 

As explained in ChapterjJl the Peres-Horodecki criterion, that is, the positivity of the partial transpose of the state, 
is a necessary condition for a state to be separable. It turns out that it is also a sufficient criterion for 1 x 1 Gaussian 
states. It is possible to show that for continuous variables, transposition is equivalent to a change of sign of the 
momenta; hence, the partial transposition operation corresponds to a change of sign in the corresponding momenta. 
In the case of a I x I Gaussian state this means that partial transposition of the second mode is equivalent 

to the transformation 

r = {l2x2®2)f, V = {12x2® Z)V{12x2 ®Z). (6.147) 

The Peres-Horodecki criterion is then reduced to check whether H is a physical covariance matrix. It is not difficult to 
prove that V is positive definite, and hence, the only condition left to analyze is H -I- iH > 0, or, equivalently, > 1 
in terms of the symplectic eigenvalues of V (note that being symmetric and real, V satishes Williamson’s theorem as 
well). 

Let us consider the example of the two-mode squeezed vacuum state |r)i 2 , for which 


V = 


1 - 2 x 2 cosh2r 
— 12 x 2 sinh2r 


— 12 x 2 sinh2r 
12 x 2 cosh 2r ’ 


which has 

det1d = l, and A(1d) = 2(cosh^ 2r-|-sinh^ 2r) = 2(1-|-2 sinh^ 2r), 
and therefore symplectic eigenvalues 


(6.148) 

(6.149) 


iy± = {l + 2 sinh^ 2r) ± ^{1 + 2 sinh^ 2r)2 - 1. (6.150) 

For any r > 0 we have z/_ < 1, what is a signature of |r)i 2 being an entangled state. 

From an experimental point of view, the Peres-Horodecki criterion is quite demanding, as it requires the full 
reconstruction of the covariance matrix of the observed beams. However, there is a simpler separability criterion which 
requires only the analysis of the variance of a suitable pair of joint quadratures (what can be check experimentally 
via two homodyne measurements, as we will see later). This criterion, which was introduced simultaneously by Duan, 
Giedke, Girac, and Zoller m and by Simon [T2], states that a Gaussian state is separable if for every ip 

= V[{X^ - Xl)l^f2\ + V[{P^ + P^)/V2] > 2, (6.151) 


where 


^ givgt _|_ g = Xcosi^ + Psin(/j, (6.152) 

= i(e“^a’^ — e“''^a) = P cos ip — X sin(/?, (6.153) 

and V{A) = {SA^). It is possible to show that this is also a sufficient condition for separability in the case of 1 x A 
Gaussian states. Note that for covariance matrices written in the standard form (16.1431) . the so-called witness 
reduces to 

= a-I-5-I-(c2 — Cl) cos2(/5. (6.154) 

Let’s come back to the example of the two-mode squeezed vacuum state |■r)l 2 ; its covariance matrix (16.1161) being 
already in standard form, we get 

= 2(cosh2r-I-sinh2r cos2(/j); (6.155) 

for ip = 7r/2, the witness reads 

= 2exp(-2r), (6.156) 

which is clearly below 2 for every r > 0, hence showing once again that |r)i 2 is indeed an entangled state. In 
particular, |r)i 2 presents quantum anti-correlations between the position of the oscillators, and correlations between 
their momenta, that is. 


^[(A - A)/A2] = H[(Xi + X 2 )IV 2 ] = exp(-2r); 


( 6 . 157 ) 
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in the r —>■ oo limit, the state |r)i 2 is therefore an eigenstate of the Pi — P 2 and Xi + X 2 operators given by 


|EPR) = j dp\p,p) = j dx\x,—x)., (6.158) 

Jr Jr 

which is exactly the type of state that Einstein, Podolsky, and Rosen considered in his attempt at proving that 
quantum mechanics was inconsistent m- 

Let us finally stress that a necessary and sufficient criterion for separability has been found for the Gaussian states 
of a general N x M bipartite continuous variable system. This criterion states that the Gaussian state pg{y^V) is 
separable if and only if there exists a pair of matrices Va and Vb with dimensions 27V x 27V and 2M x 2M, respectively, 
for which 

V>Va(SVb. (6.159) 

Of course, this criterion is quite difficult to handle in practice, but fortunately an equivalent, operationally friendly 
criterion was introduced by Giedke, Kraus, Lewenstein, and Cirac, based on the concept of nonlinear maps; we will 
however not explain this criterion which can be consulted in |14) , or in the original reference |15) . 

Let’s move now to the quantification of the entanglement present in a Gaussian state. As we commented in 
Ghapter [31 this problem has been only solved for pure states, for which the entanglement entropy is the unique 
measure of quantum correlations. For mixed states, however, we have not found a completely satisfactory measure 
(the distillable entanglement and the entanglement of formation are reasonable measures, but cannot be computed for 
most states, while the logarithmic negativity is easy to compute but does not satisfy all the conditions needed for a 
proper entanglement measure), not even for the reduced class of continuous variable Gaussian states. 

In the following we explain how to compute the entanglement entropy and the logarithmic negativity for 1x1 
Gaussian states poijy f^)i whose covariance matrix we write in the same block form as before 


V = 


A 


C 

B ’ 


(6.160) 


where A = A^, B = and C are 2x2 real matrices (the generalization to NxM Gaussian states is straightforward). 

As commented in (16.671) . tracing out one mode of a Gaussian state is equivalent to retaining in the covariance matrix 
the minor corresponding to the remaining modes. Hence, to evaluate the entanglement entropy of the Gaussian state 
Pg(j^ y) having covariance matrix (I6.16()|) . one just needs to evaluate the entropy of the single-mode covariance matrix 
A. This matrix has v = VdetA as its sole symplectic eigenvalue, and therefore, based on (16.13611 . its entropy—and 
hence the entanglement entropy of the corresponding two-mode Gaussian state—is given by 

E[pG{f, V)] = g{Vdet A), (6.161) 

where the function g{x) was defined in (16.1371) . 

For mixed states the entanglement entropy is not even an entanglement monotone, and hence, it cannot be consid¬ 
ered a proper entanglement measure for such states. One has then to consider other measures, and here we focus on the 
logarithmic negativity E]v[p]. It is possible to show that for an arbitrary Gaussian state /3 g(F, E), this entanglement 
measure can be computed as 

EN[pG{f,V)]=J2E{i>,), (6.162) 

3 

where 

= { , (6.163) 

and is the symplectic spectrum of the covariance matrix V corresponding to the partial transposition of /3G(f, V), 
which is defined in (16.1471) for a 1 x 1 system. 

Let us start by quantifying the entanglement of the two-mode squeezed vacuum state |r)i 2 ; being a pure bipartite 
state, its entanglement is measured by its entanglement entropy 


A[|r)i 2 ] = g(cosh2r), 


(6.164) 


which is nothing but the entropy of the reduced thermal state. It is not difficult to check that, starting at zero for 
r = 0, this is a monotonically increasing function of r, as expected. 

In the case of the two-mode squeezed vacuum state we can even prove a stronger result, namely that |r')i 2 1^ |r)i 2 
for r' <r. Let us write the state as 

00 

|A)i 2 = X! \/Pn(A)|n,n), 

n=0 


(6.165) 
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where 

p„(A) = (1 — being A = tanhr. 

It is fairly simple to check that the triangular matrix 


DiX,X') 


do 0 0 0 

di do 0 0 

d2 di do 0 

da d2 di do 


with d„ = j- H{n- 1)A'^]A^(” 
1 — A'^ 


(6.166) 


(6.167) 


H{x) being the Heaviside step function defined as H{x) = 1 for x > 0 and H{x) = 0 for x < 0, is column stochastic, 
and that it connects the Schmidt distributions p(A) and p(A') as 

p(A) = £i(A,A')p(A'). (6.168) 

Hence, as p(A') can be transformed into p(A) via a column stochastic matrix, we conclude that |A')i 2 majorizes 
|A)i 2 for X' < X. This implies the previous result that we found concerning the entanglement entropy, namely that 
£'[k)i 2 ] > 7i^[|»’)i2] for r > r', and much more, for example, that |r')i 2 can be transformed into a two-mode squeezed 
vacuum state of lower entanglement via an LOCC protocol. 

6.3.5 Gaussian channels 
General definition 

In this section we introduce one of the most important objects in the field of quantum information with continuous 
variables: Gaussian channels. We call channel to any trace preserving quantum operation acting on a continuous 
variables system. The channel is Gaussian when it maps Gaussian states into Gaussian states. As we will see, 
they receive their name because they actually model the most important communication channels used in current 
technologies, such as fibers or wires. 

As we saw, a way of characterizing an arbitrary trace preserving quantum operation £ is by giving a complete set 
of Kraus operators {Ek}k=i, 2 ,...,K which transform a state p into the state 

K 

p' =£[p] = Y,Ek0^. (6.169) 

fc=i 

Gaussian channels, on the other hand, can be characterized by their action on the first moments of the state. In 
particular, similarly to Gaussian unitaries, and as will be clear from the following discussion, Gaussian channels acting 
on N modes of the electromagnetic field are characterized by a vector d e plus two real 27V x 27V matrices K. 
and TV, which transform the mean vector f and covariance matrix V of the state as 

r^/Cr-Hd, and V ^ICVIC^ +M. (6.170) 

The matrices K, and TV must satisfy certain conditions in order to correspond to a true trace preserving quantum 
operation. First, as the covariance matrix is symmetric, so has to be the matrix TV. Secondly, in order to map positive 
operators (like density matrices) into positive operators, they have to satisfy the following restriction 

TViH - i/CH/C'^ > 0. (6.171) 

Note that Gaussian unitaries correspond to a Gaussian channel for which TV is zero, and K. is symplectic. Note also 
that for single-mode channels (TV = 1), this last condition can be rewritten as 

TV>0, detTV> (det/C-1)V (6.172) 

Roughly speaking, K, plays the role of the amplification and attenuation of the channel (plus a possible rotation), while 
TV includes any source of quantum or classical noise; we will come back to their physical meaning in the next section. 

Indeed, it is quite simple to understand why Gaussian channels correspond to a transformation of the type 1)6.1701) . 
To this aim, we just need to remember that any trace preserving quantum operation can be seen as a unitary 
transformation acting on the system, plus some environment in a pure state which is dismissed after the interaction. 
It is obvious that in order for the channel to be Gaussian, both the state of the environment \iPe) and the joint unitary 
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transformation Ug{s,S) must be Gaussian. Moreover, as every pure Gaussian state is connected to vacuum via some 
Gaussian unitary transformation which can be included in the joint unitary Ug(s,S), we can take the initial state of 
the environment as the multi-mode vacuum state, that is, 


Ne 

\'<Pe) = (^|0) = |uac). 




(6.173) 


where we have assumed that the environment consists in Ne modes. The state p of the system is then transformed 
into 

p' = tr£;{f/G(s,5)[/5 0 |uac)(uac|]t/^Q(s, 5)}. (6.174) 

Let us write the Gaussian parameters associated to the joint unitary as 


s = col(d, ds). 


and 



SsE 
Se ’ 


(6.175) 


where d S and the real matrices 5s, SE^ Ssei and SeSi have dimensions 27V x 27V, 2Ne x 2Ne, 27V x 2Nej and 
2Ne X 27V, respectively. Let us write also the mean and the covariance matrix of the initially separable joint state of 
the system plus the environment as 


rs_E = col(r, 0), and Vse = 

V ®I2Nex2Nb', 

(6.176) 

after the unitary, these are transformed into 



f'sE = SfsE = col(5sr -t- d, d^). 


(6.177) 

y! — — ['^•5 ^*^5 + SseSse 

Vse-SVseS 

SsVSs + SseSe 

SesVSes + SeSe 

(6.178) 


so that by tracing out the environment, the transformation onto the mean vector and the covariance matrix of the 
system is 

r' = 5sf -f d, and V' = SsVSg + SseSse, (6.179) 

which is exactly the type of transformation introduced in (16.17011 . where we now make the identifications 

IC = Ss, and J\f = SsESgE- (6.180) 


We would like to stress that it is possible to show that the Stinespring dilation of any Gaussian channel can be 
generated by choosing an environment with less than twice the number of modes of the system, that is, Ne < 27V. 

Note that the Gaussianity is a property of the channel, not of the state of the system, that is, one can consider 
the action of the Gaussian channel onto non-Gaussian states, as we shall make later. Indeed, the transformation of a 
general state p after passing through the channel receives a very simple description in terms of characteristic functions. 
To see this, let us write the inverse of the symplectic matrix S as 


'Ts 

Tse 


{S/Se)-^ 

-Ss^Sse{S/Se)-^' 

Tes 

T'e. 


_-{S/Se)-^SesSs^ 

[S/Ss)-^ 


where we have used the general block-inversion formula (see the Wikipedia!), being 

S/Ss = Se — SEsSg^SsE and S/Se = Ss — SseSe^Ses, 


(6.181) 


(6.182) 


the so-called Schur complements of Ss and Se, respectively. Let us also denoting by r and the phase space 
coordinates of the relevant modes and the environmental modes, respectively, so that the initial characteristic function 
can be written as 

Xo{r,rE) =XpiTt)Xvac(.rE); (6.183) 

then, recalling the transformation of the characteristic function under Gaussian unitaries (I6.84|) . and after tracing 
out the environmental modes (that is, setting to zero their phase space variables), we get the output characteristic 
function 

X£[p](i') = e~'^’^^^'^XpiTsr)xvaciTEsr). (6.184) 

In the following, we will denote by £{IC,Af) any Gaussian channel, obviating the displacement d which can actually 
be generated after the channel via a unitary displacement transformation (j6.90|l . since this does not change the 
covariance matrix in any way. 
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An example: phase-insensitive Gaussian channels 

There is a particularly simple class of single-mode Gaussian channels which play an important role in communication 
technologies: the phase-insensitive Gaussian channels, which are defined by 

/C = A/rl 2 x 2 , and N = pX 2 x 2 , (6.185) 


where r > 0 and p > 0 satisfy 

/i > |t — 1|, (6.186) 

by virtue of the positivity condition (I6.172|) . Note that this class of channels are called “phase-insensitive” because 
they are invariant under rotations in phase space. We will denote them by C(t,p). 

After crossing the channel, the covariance matrix V of any state is transformed into 


V' = KVK,^ +N = 


rVii-h p tVi2 
TV 21 tVii -I- p 


(6.187) 


and hence 


trV' = TtYV + 2p, (6.188) 

det V' = det V -I- p[TtrV -\- p]. (6.189) 

Taking into account that the mean number of photons is proportional to the trace of the covariance matrix—see 
(16.601) — , and the von Neumann entropy is a monotonically increasing function of its determinant—see (16.1361) and 
(16.1381) — , we conclude that r acts as an attenuation (for 0 < r < 1) or amplihcation (for t > 1) factor, while p adds 
noise (mixedness) to the state. Note that this implies that quantum mechanics does not allow to attenuate or amplify 
a signal without introducing noise (at least deterministically, that is, via trace preserving operations), what comes 
from the fact that the uncertainty principle (which in turn comes from the position and momentum commutators 
themselves) must be satisfied at all times. 

There are two interesting limiting cases: 

• For 0 < r < 1 and p = 1 — t we talk about pure-loss channels, which are a good approximation of the fibers used 
in current optical communication technologies. It is fairly simple to check that the simplest Stinespring dilation 
of such channels consists in mixing the the input mode with a single environmental mode in a beam splitter 
()6.117al) with mixing angle cos^/3 = r (see Figure). The parameter T = cos^/3 is known as the attenuation 
factor or transmissivity, and the channel is usually denoted by CiT). 

• For r > 1 and p = t — 1, we talk about quantum-limited amplifiers, which are the less noisy (deterministic) 
amplifiers that quantum mechanics allow. Again, it is simple to check that the simplest Stinespring dilation 
of these channels consists in mixing the input mode with a single environmental mode in a two-mode squeezer 
(16.1071) with squeezing parameter satisfying cosher = r (see Figure). The parameter G = cosher is known as 
the amplification factor, and the channel is usually denoted by A(G). 

It is simple to see that any phase-insensitive Gaussian channel can be seen as the concatenation of a pure-loss 
channel and a quantum-limited amplifier, that is, C(r, p) = C{T) o A{G), where 

t = TG, and p = G{l-T) + {G -1). (6.190) 

6.4 Measurements in continuous variables systems 

6.4.1 General description of measurements in phase space 

In Section 321 we learned that the most general measurement that one can perform in a quantum system can always be 
described by a complete set of trace-decreasing operations {£j}j=i, 2 ,...,j>i, each corresponding to one of the possible 
measurement outcomes. When all the trace-decreasing operations were described by a single Kraus operator, we 
talked about a POVM-based measurement, the simplest generalizations of the familiar projective measurements. In 
this section we will learn a convenient way of describing such generalized measurements for continuous variables 
systems. 

Let’s start with some useful definitions. Given the initial state of a system p, and the POVM {flj = EjEj}j=i^ 2 ,...,j, 
we will denote by 

p^=£^[p]=E,pEl ( 6 . 191 ) 
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the unnormalized state obtained after the outcome j appears; such an outcome appears with probability pj = tr{/5j}, 
and the normalized state of the system reads pj = pj^pj. Similarly, and assuming that the system is described as a 
collection of N oscillators, we define the corresponding unnormalized characteristic and Wigner functions as 

p j2N ^ . 

Xj(r) = tr{£>(r)pj} and Wj{r)= —(6.192) 

7r2n (47r) 


from which the probability of the corresponding outcome can be obtained as 


Pj = Xi(0) = 


/r2N 




and the normalized functions as 


Xj (r) = tT{D{r)pj } = pj (r) and Wj (r) = 




(47r)2Af 




= m ^Wj{r). 


(6.193) 


(6.194) 


There are many situations in which the measurement is not applied to the whole system, but only to one of the 
modes that conform it; we talk then about partial measurements. Moreover, as we will see in the next sections, the 
measurement performed onto a light beam is usually destructive, that is, the mode disappears after the measurement 
is done, so that one has to trace it out of the system. Assuming that the system has + 1 modes, and that the 
measurement is applied to the last mode, this means that the (unnormalized) state of the remaining N modes after 
the measurement will bcEl 

Pj = trN+i{{lN 0 Ej)p{lN 0 E^)} = trAr+i{(Jjv 0 (6.198) 

where p is the initial state of the TV + 1 modes (obviously the first N modes only ‘feel’ the measurement if they share 
some correlations with the measured mode). The first interesting feature of such partial, destructive measurements is 
that one only needs the POVM {tlj}j=i^ 2 ,...,j to evaluate the final state of the non-measured modes; this is contrast 
to non-destructive measurements, which require knowledge of the measurement operators {£'j}j=i, 2 ,...,j in order to 
understand the final state. 

Partial measurements iare easily described in terms of the characteristic and Wigner functions. In the case of the 
characteristic function, the derivation is simple by using both sides of (16.4911 : 

Xj{^{N}) = triN}{D{r{N})Pj} = tr{[L)(r|Ar}) 0 flj]p} (6.199) 

Xp(r)xn,-(-i~^-n)- 

(47r)^(5(2A'‘) (r^iv}—s{iv>) 


p (i2(iV+l)g ^ _ p 

= / NAf+i Xp(s)ti'{iv}{J>(r{Af} - S{jv|)}trjv+i{nj-j:)T(sjv+i)} = / 

7r 2(N- + 1) (47r) + V-^-A--/ Jr2 


Using (I6.192|) and (16.5011 . we derive now the transformation rule for the Wigner function: 

d'^^S{N} f d‘^SN+1 


Wj{r{N}) = 


xMN},SN+i)Xn, (-SAr+i)ei^{->''--t-> 


R2W 


(47r) 


2N 


= / d^TN+l 


47r 

rf2(iV+l) 


( 6 . 200 ) 


rR2(Ar + l) 


7;r^;^l^Xp(s)Wn (rAr+i)e = ® ’^ = An d rN+iWpir)Wf^(rN+i). 

^ ^ j^2 J 


Hence, both for the characteristic and the Wigner functions, the transformation is obtained by multiplying the initial 
function of the N 1 modes by the function associated to the POVM element Ilj (with the suitable sign in the 
argument), and integrating out the measured mode. 

® Note that the cyclic property applies also for the partial trace when the operator acts as the identity on the non-traced subspaces. In 
particular, consider a bipartite Hilbert space space "Ha and let us denote by {|aj)}j and orthonormal bases in the individual 

subspaces. Taking this into account, let’s prove the following identity: 

tr{/5(JA t8) S)} = tr{(/A (8)-B)p}, (6.195) 

where p is an operator acting on the complete space "Ha 0 'Hb- To this aim, we just use an explicit representation of the p operator 

p=E Pkl,Tnn \a.k){ai\^\bm){bn\ tr{/3(/A ® S)} = EE Pkl,Tnn \a.k){ai\{{bj\bm){bn\B\bj))-, (6.196) 

klmn j klmn 


now, introducing the identity Ib I^T’){^r| right before B, we get 

tr{p(/A 8) S)} = y y Pkl,mn\uk){ai\ (0 ( X) |bj) I | \bm){bn\br) 

r klmn \ j J 


— ® ^) ( E Pkl,mn \uk){ai\ ® \hm){bn\ I |br), 

(6.197) 


which is exactly the right had side of 116.19511 . 
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Fig. 6.1: Schematic representation of Mollow’s ideal single-mode detection. 


6.4.2 Photodetection: measuring the photon number 

The most fundamental measurement technique for light is photodetection. As we shall see with a couple of examples 
(homodyne detection and on/off detection), any other scheme used for measuring different properties of light makes 
use of photodetection as a part of it. 

This technique is based on the photoelectric effect or variations of it. The idea is that when the light beam that we 
want to detect impinges a metallic surface, it is able to release some of the bound electrons of the metal, which are then 
collected by an anode. The same happens if light impinges on a semiconductor surface, though in this case instead of 
becoming free, valence electrons are promoted to the conduction band. The most widely used metallic photodetectors 
are known as photo-multiplier tubes., while those based on semiconducting films are the so-called avalanehe photo 
diodes. In both cases, each photon is able to create one single electron, whose associated current would be equally 
difficult to measure by electronic means; for this reason, each photoelectron is accelerated towards a series of metallic 
plates at increasing positive voltages, releasing then more electrons which contribute to generate a measurable electric 
pulse, the photopulse. 

It is customarily said that counting photopulses is equivalent to counting photons, and hence, photodetection is 
equivalent to a measurement of the number of photons of the light field. This is a highly idealized situation, valid only 
in some limits which we will try to understand now. 

Consider the following model for a perfectly efficient detection scheme. A single-mode field with boson operators 
{a, a^} initially in some state p is kept in continuous interaction with a photodetector during a time interval T. The 
intuitive picture of such a scenario is shown in Figure [53] A cavity formed by the photodetector itself and an extra 
perfectly reflecting mirror contains a single mode. By developing a microscopic model of the detector and its interaction 
with the light mode, Mollow was able to show that the probability of generating n photoelectrons (equivalently, the 
probability of observing n photopulses) during the time interval T is given by [Iti] 


Pn = 






exp [—(1 — e '‘^)a’^a] : \ , 


( 6 . 201 ) 


where the expectation value has to be evaluated in the initial state p of the light mode, and k is some parameter 
accounting for the light-detector interaction. Using the operator identity :exp [—(1 — e~^)d)ii\: = exp(—Aa’^’d) [T7] . 
and the help of the number state basis {|n)}„gN, it is straightforward to get 




to ! 


Pn = > {n\p\n)— -- 

n\[m — n)\ 

m—n ^ ' 




-kT\ 




{n\p\n), 


( 6 . 202 ) 


and hence, for large enough detection times the number of observed pulses follows the statistics of the number of 
photons. In other words, this ideal photodetection scheme is equivalent to measuring the number operator a^a as 
already commented, that is, a projective measurement with projectors {P„ = |n)(n|}„=o.i, 2 ,...- 

However, in real photodetectors the condition T ^ is hardly met; one usually dehnes the quantum effieieney 
r] = 1 — e~'^'^, which in current photodetectors varies from one wavelength to another, and then photodetection is 
equivalent to a generalized measurement with POVM elements {n„}„=o,i,...! being 
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This POVM-based measurement admits a very simple Stinespring dilation, which is quite convenient to gain some 
intuition about optical measurement schemes: before arriving to a photodetector with unit quantum efficiency, the 
optical mode is mixed with an ancillary vacuum mode in a beam splitter of transmissivity cos^ (3 = ij. In order to 
prove that this scheme leads to the same POVM as the detector with finite efficiency, let us compute the probability 
of observing n photopulses in the detector. Using the identity 


_ ^l3aa\,-l3a'^aE _ gOot^ 


B{/3) = 


tan (3 


(cos py 


and taking into account that 




|fc,0) 


k 

E 

3=0 


'ti 


|A:,0)=^C^ 
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(6.205) 


the state of the system after the beam splitter can be written as 


PSE = B{/3){p<S)\0){0\)B\/3) = cos"+™/3^ ^ . 

n,m—0 j=0 l—O 


n\ m 


taffi+'/3|n-j,j)(m-/,Z|. 


(6.206) 


The reduced of the detected mode reads then, 

OO 

/OS = tr{ps_E} = Y, Pnm COS 
n,771=0 


min{7i,77i} 
.n+m p ^ 

fc =0 


n\ f m 


tan^'^ /3\n — k){m — A:|, 


(6.207) 


and since the detector is taken as ideal, the probability of observing I photopulses is equal to 

OO 

tan PSn — k,l^m — k,l — ^ ^ P 


OO min{?i,m} 

Pl = {l\ps\l)= ^ ^ ^ 


71,771 = 0 


fc = 0 


n \ / m 

kl[k 


SnmSk.n-l 


n—l 


n — I 


cos2"/3tan2(’"-')/3, 
(6.208) 

which coincides with (16.2021) once the identification cos^ /3 = ry is done. 

Apart from the finite quantum efficiency, which accounts for the missed photons which do not generate photoelec¬ 
trons in the detector, there is another source of imperfection in the photodetector: electrons which are pulled out from 
the detector without interacting with any photon of the detected mode. One refers to the corresponding photopulses 
as dark counts, and they can be modeled within the previous Stinespring dilation in a very simple way: by assuming 
that the ancilla mode is not in vacuum but in some other state, say pE. In this scenario, the POVM elements become 


n„ = tr{(J (g) pE)B\mn){n\ 0 i)B{3)}. 


(6.209) 


In the following all these imperfections will be ignored, so that we will assume that photodetection is equivalent to 
a measurement of the number of photons of the field impinging the detector. However, it is important to understand 
the experimental limitations before proposing any interesting theoretical protocol, and also to know how to treat them 
theoretically in case we find the need of doing a more realistic analysis. 


6.4.3 Homodyne detection: measuring the quadratures 

Even though the output of the photodetectors can take only integer values (number of recorded photopulses), they 
can be arranged to approximately measure the quadratures of light, which we remind are continuous observables. This 
arrangement is called homodyne detection. The basic scheme is shown in Figure IHIH The mode we want to measure 
is mixed in a beam splitter with another mode, called the local oscillator, which is in a coherent state |aLo)- When 
the beam splitter is 50/50 the homodyne scheme is said to be balanced, and the annihilation operators of the modes 
leaving its output ports are given by 

a± = ^(a±aLo), (6.210) 

v2 

being olo the annihilation operator of the local oscillator mode. These modes are measured with independent pho¬ 
todetectors, and then the corresponding signals are subtracted. Based on the idealized photodetection picture of the 
previous section, this scheme is analogous to a measurement of the photon number difference 

No = a+a+ — d^_d- = aj^Qa -I- OLoh^. 


(6.2II) 
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Fig. 6.2: Homodyne detection scheme with ideal photodetectors. When the local oscillator is in a strong coherent state, this 
setup gives access to the quadratures of light. 


Taking into account that the local oscillator is in a coherent state with amplitude ano = |ctLo| exp(i<p), and is not 
correlated with our measured mode, it is not difficult to show that the first moments of this operator can be written 
as 


where 


(iVD) = |aLo|(^^) 



|aLo|^ 




(a’I'd) 


= 




(6.212a) 

(6.212b) 


(6.213) 


is a generalized quadrature which coincides with the position and momentum for ip = 0 and 7r/2 respectively. Hence, 
in the strong local oscillator limit |q!loP ^ max{l, the output signal of the homodyne scheme has the mean 

of a quadrature of the analyzed mode (the one selected by the phase of the local oscillator), as well as its same 
variance. Moreover, it is simple but tedious to check that all the moments of Nd coincide with those of X‘^ in the strong 
local oscillator limit, and therefore, balanced homodyne detection can be seen as a measurement of the corresponding 
quadrature. 

It might be difficult to accept that a measurement of IVd, which has a discrete spectrum, can be equivalent to 
a measurement of X'^, which has a continuous spectrum. The reconciliation between this two pictures comes from 
the condition lapoP 1; which essentially means that the local oscillator is very intense, and therefore, there are 
so many photons impinging the detectors that the photopulses are generated at a rate much faster than the response 
time of the photodetectors, and the output signal is basically felt as a continuous photocurrent the observer. 

Just as we explained in Section 16.4.11 sometimes it is interesting to apply a measurement onto one mode out of a 
collection of modes (say the last mode of a system with iV +1 modes); this is what we defined as a partial measurement. 
Partial homodyne detection receives a very simple treatment in terms of Wigner functions. For example, in the ideal 
case explained above (strong local oscillator limit), homodyne detection of the position quadrature is described by 
the continuous set of projectors {Fix) = |a:)(a;|}a;GR; assume that the measurement pops out the outcome^ xq, so 
that, according to (16.2001) and (16.1991) . the characteristic and Wigner functions of the remaining N modes collapse to 


^ Being a measurement of a continuous observable, in real experiments the outcome cannot be a definite value; instead, one can just ensure 
that the outcome was in certain interval according to the precission of the measurement device. Nevertheless, we will use definite-outcome 
idealization, since nothing qualitatively different is introduced otherwise. 
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(unnormalized) 

/* r f 

Xxoir{N})= ^^Xp(i-)Xko)<^o|(-riv+i) and 1^^^(r{Ar}) = dvr / d^r^v+iWp(r)VF|^„)(^„l(rAr+i), 

JR2 47r 7 r2 

(6.214) 

where p is the initial state of the + 1 modes. The characteristic function of the projector |a:o)(a:o| is easily found as 


X|xo)(.o|(r) = tr{.D(r)|xo)(xo|} = e-‘P"/4(xo|e'P^/2g-ixP/2|^^^ ^ g-ipx/4gipxo/2 / dpue-'^-P'^/^xu\pu){pM = 5{x)e^P-° 


4 'k 6 {x) 


1 /47r 


(6.215) 

where we have used (|6.15l) , while Fourier transforming this expression we obtain the corresponding Wigner function 


W, 


dx'dp' 


= / 77viTX|xo>(a.o|(^'^P')e^ 

147rj 




S{x - xo) 


These expressions lead us to the following characteristic and Wigner functions of the remaining modes: 


(6.216) 


Xxo(r{ 7 V})= [ ^Xp(r{Ar},0,p)e and = f dpWp{r^f^},xo,p). (6.217) 

JR2 47r 7 r 


Note that, even though the characteristic (16.2151) and Wigner (16.2161) functions of the projector are not normalizable— 
what makes sense, since the position eigenstate isn’t either—, the characteristic and Wigner functions of the remaining 
modes can be normalized, so that the probability density function associated to the possible outcomes xq is given by 


P{xo) = XxoiO) 



d^r^N}W^g{r{M})- 


(6.218) 


Note that the projector |a:o)(xo| is Gaussian, and therefore, the partial homodyne measurement maps Gaussian 
states into Gaussian states. In the following section we evaluate the collapse of a general Gaussian state associated to 
a particular Gaussian POVM-element, and we will give there the transformation rules for such a particular class of 
states under a partial homodyne measurement. 


6.4.4 Partial on/off detection: de-Gaussification by vacuum removal 

Despite the incredible advances in photodetector technologies, practical photon-counters are still out of reach. A less 
demanding detection strategy is the so-called on/off detection, in which the detector gives a signal whenever one or 
more photons reach it, but this signal is identical no matter how many photons triggered it. In the next section we 
will see that on/off detection allows for the implementation of interesting operations, such as photon addition and 
subtraction onto a light field, which are the most fundamental non-Gaussian and non-unitary operations that one can 
think of. 

On/off detection is then a projective measurement with two possible outcomes off (‘no-click’) and on (‘click!’), 
with corresponding projectors 


A// = |0)(0| and P„n = /-|0)(0| = ^|n)(n|. (6.219) 

n—1 

These projectors have very simple characteristic functions 

Xp^^/r) = X|o)(o|(i’) and yp^^Jr) = 47r(5(2)(r) _ (6.220) 

and very simple Wigner functions as well 

W = ^loXolW and bFp„„(i-) = ^ “ W^|o)<o|(r)- (6.221) 

Applied as a partial measurement (as usual, the measurement is applied onto the last mode of a system with N +1 
modes), on/off detection makes the system evolve from some initial state p of the A^-|- 1 oscillators, to a reduced state 
of the first N oscillators with (unnormalized) characteristic and Wigner functions 


(* d^ (* 

Xoff{r{N})= ^^Xp(r)X|o)<o|(-riV-si) ^ ^o//(r{Ar}) = 47r / d2rAf+iWp(r)W|o)(o|(rAf+i) (6.222) 

JR2 47r 
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or [pn = trAr+i{p} is the reduced initial state of the N first oscillators) 

Xon{Y{N})=XpNi.^{N})-Xoff{r{N}) ^ Won{Y {N}) = Wp^{Y)-Wof f{v {N}) , (6.223) 


depending on the result of the measurement. 

Note that Wp (r) is Gaussian, and therefore, the off event projects the state of the non-measured modes into 
another Gaussian state. This is not the case for the on event, whose associated Wigner function Wp (r) is not 
Gaussian, and therefore, it can be used as a de-Gaussifying operation. In order to understand this better, let’s analyze 
the case in which the initial state is a general Gaussian state PG(d, V) of the form (16.661) with M = N and M' = 1, 
whose mean vector and covariance matrix we write as 


!■ = (i’{Ar},rAr+i), and 


V = 


V{N} C 
Vm+1 


(6.224) 


where r{Ar} S r^v+i € V{Ar} and Vjv+i are real, symmetric matrices of dimensions 2N x 2N and 2x2, 

respectively, while C is a real 2N x 2 matrix. Using the Gaussian integral (j6.55p . it is straightforward to prove that 
the probability of the off event is 


2 

Poff = /, 

y det(yAr+i + 12x2) 


1 -T 


(pAf+l + 2i2x2)rAr+i , 


(6.225) 


while the corresponding output state is the Gaussian po// = pG(ro//, Vo//) with 

To// = r{jv} — (^(VAr+i +l2x2)~^FAr+l and Voff = V{iq} + C{Vn+1 ■ (6.226) 

The probability of the on event is then pon = 1 — Poff, which, based on (16.2231) . has the associated output Wigner 
function 

Wou{V{Ny) = (l-Po//)”MW^PG(d{«}.y{iv})('^{Af})-Po//'^PG(d„//.y„//)(l’{Ar})] (6.227) 

Note that even though this Wigner function is not Gaussian, it is a simple combination of two Gaussians (a “negative” 
mixture in particular), and hence this way of de-Gaussification is very convenient from the theoretical point of view, 
since all the tools of Gaussian states and operations can be used. 

Let us consider a simple example: we have two modes in the two-mode squeezed vacuum state (I6.115L and we 
perform an on/off detection onto the second mode. Given the photon number correlation between the modes, it is 
obvious that whenever the outcome is off, the first mode gets projected into the vacuum state, poff = |0)(0|; on the 
other hand, if the outcome is on the state of the first mode will collapse to the (unnormalized) mixture 


OO 

Ron = y^tanh"r|n)(n|. 

n=l 


(6.228) 


According to (I6.225|) . and using the mean vector and covariance matrix of the two-mode squeezed vacuum state ()6.116D . 
the probability for the off event is Poff = 1/cosher, and it is simple to check that (16.2271) leads to Toff = 0 and 
Voff =1.2x2, corresponding to the vacuum state. The probability of the on event reads then pon = tanh^ r, and the 
corresponding Wigner function is 

Woniri) = sinh"^r[coshVW^^^(,i„h 2 r)(i'i) - IU|o)(o|(ri)], (6.229) 


that is, a “negative mixture” of a thermal and a vacuum state. Since the weight of the vacuum state is always larger 
than that of the thermal state, this Wigner functions always has a negative central region, surrounded by a positive 
one, what shows the non-Gaussian character of the state. Moreover, this central negative region has more or less the 
same size irrespective of the squeezing value; indeed, it is very simple to show that the radius of the central negative 
region is 


i?_ = 


'21ogg(l -l-tanh^ r) 
tanh^ r 


(6.230) 


which is a monotonically increasing function of the squeezing, but is lower and upper bounded by 1 and log^^ 4 Ri 1.18, 
and hence varies very little with r. In contrast, the positive region gets larger as the squeezing increases, what comes 
from the thermal component of the state. Finally, note that for small squeezing parameter, the state tends to the |1) 
Fock state. In fact, the next step of these notes will consist in showing that if a mode interacts very weakly via the 
two-mode squeezing or beam splitter interactions with a vacuum mode, the on detection of this second mode signals, 
respectively, the approximate application of the or d operators onto the principal mode. 

























APPENDIX 



A. THE MATHEMATICAL LANGUAGE OE QUANTUM MECHANICS: HILBERT SPACES 


Just as classical mechanics is formulated in terms of the mathematical language of differential calculus and its ex¬ 
tensions, quantum mechanics takes linear algebra (and Hilbert spaces in particular) as its fundamental grammar. In 
this section F introduce the concept of Hilbert space, and discuss the properties of some operators which will play 
important roles in the formalism of quantum mechanics. 

A.l Finite-dimensional Hilbert spaces 

In essence, a Hilbert space is a complex vector space in which an inner product is defined. Let us define first these 
terms as are used in this notes. 

A complex vector space is a set V, whose elements will be called vectors or kets and will be denoted by {|a), |6), |c),...} 
(a, b, and c may correspond to any suitable label), in which the following two operations are defined: the vector addi¬ 
tion, which takes two vectors |a) and |6) and creates a new vector inside V denoted by |a) -|- |6); and the multiplication 
by a scalar, which takes a complex number a G C (in this section Greek letters will represent complex numbers) and 
a vector |a) to generate a new vector in V denoted by a |a). 

The following additional properties must be satisfied: 

1. The vector addition is commutative and associative, that is, |a) -I- \b) = \b) |a) and (|a) -|- \b)) |c) = |a) -|- 

m + \c)). 

2. There exists a null vector \null) such that |a) -I- \null) = |a). 

3. a (|a) -I- \b)) = a |a) -I- a \b). 

4. (a -I- /3) |a) = a\a) P |a). 

5. (a/3) |a) = a(/3|a)). 

6. I |a) = |a). 

From these properties it can be proved that the null vector is unique, and can be built from any vector |a) as 0 |a); 
hence, in the following we denote it simply by \null) = 0. It can also be proved that any vector |a) has a unique 
antivector |—a) such that |a) -I- |—a) = 0, which is given by (—1) |a) or simply — |a). 

An inner product is an additional operation defined in the complex vector space V, which takes two vectors |a) 
and \b) and associates them a complex number. It will be denoted by (a|6) or sometimes also by (|a), \b)), and must 
satisfy the following properties: 

1. (a|a) > 0 if |a) ^ 0. 

2. {a\b) = {b\a)*. 

3. (|a), a |6)) = a{a\b). 

4. (|a), \b) -G |c)) = {a\b) + (a|c). 

The following additional properties can be proved from these ones: 

• {null \null) = 0. 

• (a|a)>IU) =a*{a\b). 

• (|a) + \b ), |c)) = (a|c) (6|c). 

• |(«IUI^ < («I«)(^IU 
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Note that for any vector |a), one can define the object (a| = (|a) ,•), which will be called a dual vector or a bra, 
and which takes a vector \b) to generate the complex number (|a), |6)) G C. It can be proved that the set formed by 
all the dual vectors corresponding to the elements in V is also a vector space, which will be called the dual space and 
will be denoted by V'*'. Within this picture, the inner product can be seen as an operation which takes a bra (a| and 
a ket |6) to generate the complex number {a\b), a bracket. This whole bra-c-ket notation is due to Dirac. 

In the following we assume that any time a bra (a| is applied to a ket |6), the complex number {a\b) is formed, so 
that objects like |fe)(a| generate kets when applied to kets from the left, (|6)(a|) |c) = ((a|c)) |6), and bras when applied 
to bras from the right, (c| (|fc)(a|) = ((c|6)) (a|. Technically, |6)(a| is called an outer product. 

A vector space equipped with an inner product is called an Euclidean space. In the following we give some important 
definitions and properties which are needed in order to understand the concept of Hilbert space: 

• The vectors {|ai), [02),..., \am)} are said to be linearly independent if the relation ai |ai)+Q:2 |a2) + ...+Q:m \cim) = 

0 is satisfied only for ai = 02 = = otm = 0, as otherwise one of them can be written as a linear combination 

of the rest. 

• The dimension of the vector space is defined as the maximum number of linearly independent vectors, and can 
be finite or infinite. 

• If the dimension of an Euclidean space is d < 00 , it is always possible to build a set of d orthonormal vectors 

^ 2 d satisfying {ej\ei) = Sji, such that any other vector |a) can be written as a linear superposition 

of them, that is, |a) = being the a^’s some complex numbers. This set is called an orthonormal 

basis of the Euclidean space V, and the coefficients Oj of the expansion can be found as Oj = (e^ja). The column 
formed with the expansion coefficients, which is denoted by col (oi, 02 ,..., o^), is called a representation of the 
vector I a) in the basis E. 

Note that the set E+ = ^ d orthonormal basis in the dual space V"*", so that any bra (a| can be 

expanded then as (a| = (®il- "The representation of the bra (a| in the basis E corresponds to the row 

formed by its expansion coefhcients, and is denoted by (a*, a^,..., a*). Note that if the representation of |a) is 
seen as a d x 1 matrix, the representation of (a| can be obtained as its 1 x d conjugate-transpose matrix. 

Note finally that the inner product of two vectors |a) and \b) reads {a\b) = when represented in the 

same basis, which is the matrix product of the representations of (a| and |5). 

For finite dimension, an Euclidean space is a Hilbert space. However, in most applications of quantum mechanics 
(and certainly in quantum optics), one has to deal with infinite-dimensional vector spaces. We will treat them after 
the following section. 


A.2 Linear operators in Unite-dimensional Hilbert spaces 

We now discuss the concept of linear operator, as well as analyze the properties of some important classes of operators. 
Only finite-dimensional Hilbert spaces are considered in this section, we will generalize the discussion to infinite¬ 
dimensional Hilbert spaces in the next section. 

We are interested in maps L (operators will be denoted with throughout the notes) which associate to any 
vector \a) of a Hilbert space TL another vector denoted by L |a) in the same Hilbert space. If the map satisfies 

L{a\a)+d\b)) = aL\a)+dL\b), (A.l) 

then it is called a linear operator. For our purposes this is the only class of interesting operators, and hence we will 
simply call them operators in the following. 

Before discussing the properties of some important classes of operators, we need some definitions: 

• Given an orthonormal basis E = {\ej)} ^ d ^ Hilbert space TL with dimension d < 00 , any operator L has a 
representation; while bras and kets are represented by d x 1 and 1 x d matrices (rows and columns), respectively, 
an operator L is represented by a d x d matrix with elements Lji = (le^) ,L |e;)) = (e^j L \ei). An operator L 
can then be expanded in terms of the basis E as L = i=i kt) {^i\- It follows that the representation of 
the vector \b) = L |a) is just the matrix multiplication of the representation of L by the representation of |a), 
that is, bj = J2'i=i Ljio-i- 



A.2. Linear operators in Unite-dimensional Hilbert spaces 
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• The addition and product of two operators L and K, denoted by L + K and LK, respectively, are defined by their 
action onto any vector |a): (L + K) |a) = L |a) +i^ |a) and LK |a) = L{K |o)). It follows that the representation 
of the addition and the product are, respectively, the sum and the multiplication of the corresponding matrices, 
that is, (L + K)ji = Lji + Kji and {LK)ji = LjkKki- 

• Note that while the addition is commutative, the product is not in general. This leads us to the notion of 
commutator^ defined for two operators L and K as [L,K] = LK — KL. When [L^K] = 0, we say that the 
operators commute. 

• Given an operator L, its trace is defined as the sum of the diagonal elements of its matrix representation, that 

is, tr{L} = It may seem that this definition is basis-dependent, as in general the elements Ljj are 

different in different bases. However, we will see that the trace is invariant under any change of basis. 

The trace has two important properties. It is linear and cyclic, that is, given two operators L and K, tT{L+K} = 
tr{L} + tr{if} and ti:{LK} = tr{KL}, as is trivially proved. 

• We say that a vector |/) is an eigenvector of an operator LiiL\l) = X |^); A G C is called its associated eigenvalue. 
The set of all the eigenvalues of an operator is called its spectrum. 


We can pass now to describe some classes of operators which play important roles in quantum mechanics. 

The identity operator. The identity operator, denoted by I, is defined as the operator which maps any vector 
onto itself. Its representation in any basis is then Iji = Sji, so that it can expanded as 

d 

I = (A.2) 

This expression is known as the completeness relation of the basis E; alternatively, it is said that the set E forms a 
resolution of the identity. 

Note that the expansion of a vector |a) and its dual (a| in the basis E is obtained just by application of the 
completeness relation from the left and the right, respectively. Similarly, the expansion of an operator L is obtained 
by application of the completeness relation both from the right and the left at the same time. 

The inverse of an operator. The inverse of an operator L, denoted by L~^, is defined as that satisfying 

L-^L = LL-^ = i. 

An operator function. Consider a real function / {x) which can be expanded in powers of a; as / (a;) = 
given an operator L, we define the operator function f{L) = X]m=owhere L"* means the 
product of L with itself m times. 

The adjoint of an operator. Given an operator L, we define its adjoint, and denote it by K, as that satisfying 
{\a),L\b)) = (LI |a), |6)) for any two vectors |a) and \b). Note that the representation of U corresponds to the 
conjugate transpose of the matrix representing L, that is {K)ji = L*y Note also that the adjoint of a product of two 
operators K and L is given by {KLY = UK"^. 

Self-adjoint operators. We say that H is a self-adjoint if it coincides with its adjoint, that is, H = W. A 
property of major importance for the construction of the laws of quantum mechanics is that the spectrum ^ d 

of a self-adjoint operator is real. Moreover, its associated eigenvector^ {\hj)}j=i 2 a an orthonormal basis of 
the Hilbert space. 

The representation of any operator function f{H) in the eigenbasis of H is then 
follows 

d 

f{H)=Y,f{h,)\h,){h,\. (A.3) 

^ We will assume that the spectrum of any operator is non-degenerate, that is, only one eigenvector corresponds to a given eigenvalue, 
as all the operators that appear in this thesis have this property. 


f{H) = f (hj) 6ji, from which 
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This result is known as the spectral theorem. 

Unitary operators. We say that 17 is a unitary operator if W = U~^. The interest of this class of operators is 
that they preserve inner products, that is, for any two vectors |a) and |6) the inner product (U |a), U |6)) coincides with 
(a|6). Moreover, it is possible to show that given two orthonormal bases E = {\ej)'\^ d ~ 

there exists a unique unitary matrix U which connects them as {\e'j) = U \ej)'\j=i^ 2 ,..,d, and then any basis of the 
Hilbert space is unique up to a unitary transformation. 

We can now prove that the trace of an operator is basis-independent. Let us denote by tr{L}£; the trace of an 
operator L in the basis E; the trace of this operator in the transformed basis can be written then as tT{L}E' = 
tT{U^LU}E, or using the cyclic property of the trace and the unitarity of U, tr{ii} e' = tr{[717tL} = tr{L}E, which 
proves that the trace is equal in both bases. 

Note finally that a unitary operator U can always be written as the exponential of i-times a self-adjoint operator 
H, that is, U = exp{iH). 

Projection operators. In general, any self-adjoint operator P satisfying = I is called a projector. We are 
interested only in those projectors which can be written as the outer product of a vector |a) with itself, that is. 
Pa = |a) (a|; when applied to a vector |5), this gets projected along the ‘direction’ of |a) as Pa |1>) = {{a\b)) |a). 

Note that given an orthonormal basis E, we can use the projectors Pj = \ej) {ej\ to extract the components 
of a vector |c) as Pj |c) = Cj \ej). Note also that the completeness and orthonormality of the basis E implies that 
EU Pd = ^ and P,Pi = respectively. 


Density operators. A self-adjoint operator p is called a density operator if it is positive semidefinite, that is 
(a| p |a) > 0 for any vector |a), and has unit trace. 

The interesting property of density operators is that they ‘hide’ probability distributions in the diagonal of its 
representation. To see this just note that given an orthonormal basis E, the self-adjointness and positivity of p ensure 
that all its diagonal elements {pjj}j^i 2 d either positive or zero, that is, pjj > 0 Vj, while the unit trace makes 

them satisfy Ej=i Pjj = 1- Hence, the diagonal elements of a density operator have all the properties required by a 
probability distribution. 

It is possible to show that a density operator can always be expressed as a statistical or convex mixture of projection 
operators, that is, p = l®fe) Wk\, where tUk = ^ and the vectors {|afc)}fli are normalized to one, but 

don’t need not to be orthogonal (note that in fact M doesn’t need to be equal to d). Hence, another way of specifying 
a density matrix is by a set of normalized vectors together with some statistical rule for mixing them. When only one 
vector |a) contributes to the mixture, p = |a)(a| is completely specified by just this single vector, and we say that the 
density operator is pure\ otherwise, we say that it is mixed. 

A.3 Generalization to infinite dimensions 

As explain in Chapter [51 the natural Euclidean space for quantum optical systems is infinite-dimensional, as each 
mode of the electromagnetic field behaves as a harmonic oscillator. Unfortunately, not all the previous concepts and 
objects that we have introduced for the finite-dimensional case are trivially generalized to infinite dimensions; in this 
section we discuss this generalization. 

The first problem that we meet when dealing with infinite-dimensional Euclidean spaces is that the existence of 
a basis {\ej)}j^.^ ^ which any other vector can be represented as |a) = E^=i ki) is not granted. The class of 
infinite-dimensional Euclidean spaces in which these infinite but countable bases exist are called Hilbert spaces, and 
are the ones that will be appearing in quantum mechanics. 

The conditions which ensure that an infinite-dimensional Euclidean space is indeed a Hilbert space can be found 
in, for example, reference m- Here we just want to stress that, quite intuitively, any infinite-dimensional Hilbert 
spac^l is isomorphic to the space called (oo), which is formed by the column vectors |a) = col(ai, a 2 ,...) where the 

^ An example of infinite-dimensional complex Hilbert space consists in the vector space formed by the complex functions of real variable, 
say I/) = / (x) with x G M, with integrable square, that is 

f dx\f (x) p < 00 , (A.4) 

Jr 

with the inner product 


This Hilbert space is known as the (x) space. 


(gif) = [ dxg* {x)f{x). 

Jr 


(A.5) 
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set {aj G C}j=i_ 2 ,... satisfies the restriction has the operations |a) + \b) = col(ai + 6 i, 02 + 62 ,...); 

a\a) = col(aai, aa 2 , ...)j and {a\b) = 

Most of the previous definitions are directly generalized to Hilbert spaces by taking d —i 00 (dual space, represen¬ 
tations, operators,...). However, there is one crucial property of self-adjoint operators which doesn’t hold in this case: 
its eigenvectors may not form an orthonormal basis of the Hilbert space. The remainder of this section is devoted to 
deal with this problem. 

Just as in finite dimension, given an infinite-dimensional Hilbert space 'H, we say that one of its vectors \d) is an 
eigenvector of the self-adjoint operator H if H\d) = 6\d), where J G C is called its associated eigenvalue. Nevertheless, 
it can happen in infinite-dimensional spaces that some vector |c) not contained in 'H also satisfies the condition 
H\c) = x|c), in which case we call it a generalized eigenvector^ being y its generalized eigenvalu^. The set of all the 
eigenvalues of the self-adjoint operator is called its discrete (or point) spectrum and is a countable set, while the set 
of all its generalized eigenvalues is called its continuous spectrum and is uncountable, that is, forms a continuous set 
[15] (see also [T5]l. 

In quantum optics one finds two extreme cases: either the observable, say H, has a pure discrete spectrum 
{hj}j=i^ 2 ,...', or the observable, say X, has a pure continuous spectrum {a;} 2 ;gR. It can be shown that in the first case 
the eigenvectors of the observable form an orthonormal basis of the Hilbert space, so that we can build a resolution of 
the identity as / = \^j) ^.nd proceed along the lines of the previous sections. 

In the second case, the set of generalized eigenvectors cannot form a basis of the Hilbert space in the strict sense, 
as they do not form a countable set and do not even belong to the Hilbert space. Fortunately, there are still ways to 
treat the generalized eigenvectors of X ‘as if’ they were a basis of the Hilbert space. The idea was introduced by Dirac 
[20] . who realized that normalizing the generalized eigenvectors afQ {x\y) = 5{x — y), one can define the following 
integral operator 

/ dx\x) {x\ = Jc, (A.9) 

Jr 

which acts as the identity onto the generalized eigenvectors, that is, /c|a:) = |a;); it is then assumed that coincides 
with the identity in "H, so that any other vector |a) or operator L in the Hilbert space can be expanded as 


|a) = / dxa{x)\x) and dxdyL {x,y) \x) {y\ (A.10) 

Jr Jr2 

where the elements a{x) = (a;|a) and L{x,y) = {x\L\y) of these continuous representations form complex functions 
defined in R and respectively. From now on, we will call continuous basis to the set 

Dirac introduced this continuous representations as a ‘limit to the continuum’ of the countable case; even though 
this approach was very intuitive, it lacked of mathematical rigor. Some decades after Dirac’s proposal, Gel’fand showed 
how to generalize the concept of Hilbert space to include these generalized representations in full mathematical rigor 
[2T] . The generalized spaces are called rigged Hilbert spaces (in which the algebra of Hilbert spaces joins forces with the 
theory of continuous probability distributions), and working on them it is possible to show that given any self-adjoint 
operator, one can use its eigenvectors and generalized eigenvectors to expand any vector of the Hilbert space. 

Note finally that given two vectors |a) and \b) of the Hilbert space, and a continuous basis we can use 

their generalized representations to write their inner product as 


(a|fi) 


= f dxa* (x) b (x). 

Jr 


(A.ll) 


It is also easily proved that the trace of any operator L can be evaluated from its continuous representation on 
as 



dxL [x, x). 


(A.12) 


^ In (x) we have two simple examples of self-adjoint operators with eigenvectors not contained in (x): the so-called X and P, 
which, given an arbitrary vector |/} = / (x), act as X\f) = x/ (x) and P|/} = —idf/dx, respectively. This is simple to see, as the equations 

(a;) = Xfx ( 3 ;) and - i^/p = P/p , (A.6) 

ax 

have 

fx {x) = S{x- X) and fp (x) = exp (iPx), (A.7) 

as solutions, which are not square-integrable, and hence do not belong to (x). 

^ This 5 (x) function the so-called Dirac-delta distribution which is defined by the conditions 




1 if y E [xi,X2] 
0 if y ^ [xi,X2] 


(A.8) 
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A. The mathematical language of quantum mechanics: Hilbert spaces 


This has important consequences for the properties of density operators, say p for the discussion which follows. We 
explained at the end of the last section that when represented on an orthonormal basis of the Hilbert space, its 
diagonal elements (which are real owed to its self-adjointness) can be seen as a probability distribution, because they 
satisfy Pn — ^ Pjj — 0 J- Similarly, because of its unit trace and positivity, the diagonal elements of its 
continuous representation satisfy dxp (x, x) = 1 and p{x,x) > 0 V x, and hence, the real function p (x, x) can be 
seen as a probability density function. 


A.4 Composite Hilbert spaces 

In many moments of these notes, we find the need to associate a Hilbert space to a composite system, the Hilbert 
spaces of whose parts we now. In this section we show how to build a Hilbert space H starting from a set of Hilbert 
spaces {TLa,'Hb,'Hc---}- 

Let us start with only two Hilbert spaces TIa and TIb with dimensions dA and respectively (which might be 
infinite); the generalization to an arbitrary number of Hilbert spaces is straightforward. Consider a vector space V 
with dimension dim(V) = x We define a map called the tensor product which associates to any pair of vectors 
|a) G Ha and \b) G Hb a vector in V which we denote by |a) 0 |6) G V. This tensor product must satisfy the following 
properties: 

1 . (|a) + |6)) 0 |c) = |a) 0 |c) + \b) 0 |c). 

2 . |a) 0 (|6) + |c)) = |a) 0 \b) + |a) 0 |c). 

3 . (ala)) 0 | 5 ) = |a) 0 (a|6)). 

If we endorse the vector space V with the inner product (|a) 0 |6), |c) 0 |d)) = (a|c)(6|d), it is easy to show it becomes 
a Hilbert space, which in the following will be denoted by "H = Ha 0 Hb- Given the bases Ea = {|e^)}j=i,2,...,dA 
and Eb = {|e^)}j=i,2,...,dj3 of the Hilbert spaces Ha and Hb, respectively, a basis of the tensor product Hilbert space 

Ha 0 Hb can be built as E = Ea 0 Eb = {\e-f) 0 |ef )}/=^,’2,’...,’dB (note that the notation in the first equality is 
symbolical). 

We will use a more economic notation for the tensor product, namely |a) 0 |6) = |a,6), except when the explicit 
tensor product symbol is needed for any reason. With this notation the basis of the tensor product Hilbert space is 
written as E = {\ef, 

The tensor product also maps operators acting on Ha and Hb fo operators acting on H. Given two operators 
La and Lb acting on Ha and Hb, the tensor product operator L = La 0 Lb is defined in H as that satisfying 
L|a, b) = (LaIo)) 0 (Lb|6 )) for any pair of vectors |a) G Ha and \b) G Hb- When explicit subindices making reference 
to the Hilbert space on which operators act on are used, so that there is no room for confusion, we will use the shorter 
notations La 0 Lb = LaLb, La 0 / = La, and / 0 Lb = Lb- 

Note that the tensor product preserves the properties of the operators; for example, given two self-adjoint operators 
Ha and Hb, unitary operators Ua and Ub, or density operators pA and pB, the operators Ha 0 Hb, Ua 0 Ub, and 
Pa 0 Pb are self-adjoint, unitary, and a density operator in H, respectively. Note that this doesn’t mean that any 
self-adjoint, unitary, or density operator acting on H can be written in a simple tensor product form La ® Lb- 



BIBLIOGRAPHY 


[1] M. A. Nielsen and 1. L. Chuang, Quantum information and quantum computation. Cambridge University Press, 

2000 . 

[2] R. Horodecki, P. Horodecki, M. Horodecki, and K. Horodecki, “Quantum entanglement,” Rev. Mod. Phys., vol. 81, 
pp. 865-942, 2009. 

[3] J. Eisert, “Entanglement in quantum information theory,” PhD dissertation; arXiv:quant-ph/061025. 

[41 S. L. Braunstein and P. van Loock, “Quantum information with continuous variables,” Rev. Mod. Phys., vol. 77, 
pp. 513-577, 2005. 

[5] C. Weedbrook, S. Pirandola, R. Garcia-Patron, N. J. Cerf, T. C. Ralph, J. H. Shapiro, and S. Lloyd, “Gaussian 
quantum information,” Rev. Mod. Phys., vol. 84, pp. 621-669, May 2012. 

[6] C. Navarrete-Benlloch, Contributions to the quantum opties of multi-mode optieal parametric oscillators. PhD 
dissertation, 2011. 

[7] B. Arnold, Majorization and the Lorenz order. Springer-Verlag Lecture Notes in Statistics 43, 1987. 

[8] W. P. Schleich, Quantum opties in phase spaee. Wiley-VCH, 2001. 

[9] R. Simon, N. Mukunda, and B. Dutta, “Quantum-noise matrix for multimode systems: U(?r) invariance, squeez¬ 
ing, and normal forms,” Phys. Rev. A, vol. 49, pp. 1567-1583, Mar 1994. 

[10] C. C. Gerry and P. L. Knight, Introduetory quantum optics. Cambridge University Press, 2005. 

[11] L.-M. Duan, G. Giedke, J. 1. Cirac, and P. Zoller, “Inseparability criterion for continuous variable systems,” Phys. 
Rev. Lett, vol. 84, p. 2722, 2000. 

[12] R. Simon, “Peres-horodecki separability criterion for continuous variable systems,” Phys. Rev. Lett., vol. 84, p. 
2726, 2000. 

[13] A. Einstein, B. Podolsky, and N. Rosen, “Gan quantum-mechanical description of physical reality be considered 
complete?” Phys. Rev., vol. 47, pp. 777-780, 1935. 

[14] A. Ferraro, S. Olivares, and M. G. A. Paris, “Gaussian states in continuous variable quantum information,” 
lecture notes; arXiv:quant-ph/0503237. 

[15] G. Giedke, B. Kraus, M. Lewenstein, and J. 1. Cirac, “Entanglement criteria for all bipartite gaussian states,” 
Phys. Rev. Lett, vol. 87, p. 167904, 2001. 

[16] B. R. Mollow, “Quantum theory of field attenuation,” Phys. Rev., vol. 168, pp. 1896-1919, 1968. 

[17] W. H. Louisell, Quantum statistieal properties of radiation. John Wiley & Sons, 1973. 

[18] E. Prugovecky, Quantum mechanics in Hilbert space. Academic Press, 1971. 

[19] A. Galindo and P. Pascual, Quantum Mechanics 1. Springer Verlag, 1990. 

[20] P. A. M. Dirac, The principles of quantum mechanics. Oxford university press, 1930. 

[21] 1. M. Gelfand and N. Y. Vilenkin, Generalized Funetions, Vol. IV. Academic Press, 1964. 


